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PASSIVE  INFRARED  SURVEILLANCE:  NEW  METHODS  OF  ANALYSIS 


I .  INTRODUCTION 

A.  Background 

This  paper  describes  an  analytical  model  for  predicting 
the  performance  of  a  particular  class  of  Infrared  sensors, 
generlcally  described  as  Infrared  Search  and  Truck  (IRST) 
devices.  An  IRST  system  generally  con*  ..sts  of  one  or  more 
photodetectors  located  In  the  focal  plane  of  a  scanning 
optical  telescope,  and  a  complement  of  signal  processing 
electronics  to  process  the  detected  photocurrents.  The 
signal  processors'  task  is  to  determine  whether  or  not  an 
object  of  a  particular  type  (a  "target")  Is  anywhere  in  the 
sensors'  field  of  view,  while  keeping  the  frequency  of  false 
target  reports  (l.o.,  the  "False  Alarm  Rate,"  PAR)  to  an 
acceptably  low  level. 

The  IRST  is  a  nonimaging  device,  as  contrasted  with 
Forward  Looking  Infrared  (FLIR)  imaging  systems.  ^ 

The  search  and  tr.ack  device  may  be  required  to  keep  a  full 
hemisphere  (Bti  steradians)  under  constant  observation,  have 
a  resolution  of  one  mllliradian  or  less,  and  operate  without 
human  assistance  in  a  fully  autonomous  mode  for  weeks  at  a 
time.  The  challenge  this  presents  to  the  system  designer 
is  magnified  still  further  by  the  abundant  opportunities 
for  target/background  confusion  offered  by  such  typical 
background  scenes  us  cloudy  skies  and  cities. 

VUnuMript  tubmitted  July  A,  1079, 
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Previous  attempts  to  model  background  effects  on  IRST 
system  performance  have  focused  on  the  Wiener  spectrum 
approach,  frequency  domain  technique  originally 

developed  for  calculating  the  noise  variance  in  communication 
circuits.  Unfortunately,  unlike  the  noise  processes  typically 
assumed  in  statistical  communication  theory^  the  IRST 
photocurrent  is  apt  to  be  a  highly  non-stationary  random 
process,  inadequately  characterized  simply  in  terms  of  its 
variance.  As  the  sensor  scans  across  a  structured  background, 
the  spatially  non-uniform  scene  brightness  is  mapped  Into  a 
photocurrent  whose  mean,  variance,  and  other  statistics, 
are  all  functions  of  time.  The  inadequacy  of  the  Wiener 
spectrum  method  under  these  conditions  has  been  appreciated 
for  many  years . ^ 

The  approach  taken  in  this  paper  brings  together  and 
generalizes  some  results  from  the  theories  of  optical 
communications^  “^and  stochastic  processes.^  *  In  particular, 
the  formulation  for  "threshold  crossing  rates"  originally 
derived  by  Ric^^^^for  stationary  processes  ia  generalized 
in  Section  III.B  to  accomodate  such  nonstationary  processes 
as  arise  when  nonuniform  scenes  are  scanned  by  an  IRST 

sensor,  A  crossing  rate  formulation  similar  to  that  of 

\ 

Section  III.B  has  been  previously  derived  by  Cramer  and 
Leadbetter;  '  as  discussed  in  Appendix  D. 
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The  general  croHsing  rate  formalism  is  particularized 
to  non-stationary  Gaussian  statistics  In  Section  III.C, 
appropriate  to  the  description  of  the  photocurrents  of 
present  interest! 

Some  asymptotic  approximations  to  the  complete  crossing 
rate  integral  are  developed  in  Section  III.D. 

Numerical  examples  are  given  in  Section  V  that 
compare  a  number  of  approximate  formulations  with  the 
rigorous  formulation  developed  in  Section  III.C.  An 
original  analysis  of  adaptive  threshold  IRST  sensors  is 
presented  in  Section  VI. 

Previous  works  describing  analysis  techniques  for 
IRST  systems  are  restricted  in  applicability  to  uniform 
scenes.  Moreover,  these  works  generally  prescribe  a 
totally  different  formulation  for  the  calculation  of 
FAR  than  tor  the  calculation  of  the  target  detection 
probability,  original  contribution,  of  this 

paper  is  two-told: 

1.  An  original  analysis  is  presented  that  unifies 
the  previous  work,  by  showing  that  the  earlier 
techniques  for  calculating  FAR  and  Pj^  are  actually 
both  derivable  as  special  cases  of  a  more  general 
foi'mulation ,  eq  .  (51). 

2.  More  importantly,  an  original  method  for  calculating 
FAR  and  Pp  is  presented  that  is  inherently  applicable 

to  non-uniform  ;.ceueB. 
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A  Basic  Threshold  Receiver 


In  order  to  provide  a  frame  of  reference  for  the 
following  discussion  it  is  necessary  to  describe  a  simple 
IRST  receiver  structure,  and  to  define  the  parameters  used 
to  characterize  IRST  performance. 

A  basic  threshold  oomparison  receiver  is  shown  in 
Fig.  1.  The  current  X(t)  at  the  output  of  the  detector  is 
input  to  an  electrical  filter  of  transfer  function  H(f). 

The  output  current  Y(t)  of  the  electrical  filter  is  compared 
with  a  threshold  y^Ct).  If  Y(t)  exceeds  the  threshold,  the 
presence  of  a  "target"  is  declared;  otherwise,  no  target 
declaration  is  made, 

The  probability  that  the  IRST  device  makes  a  target 
declaration  when  a  target  is  in  fact  in  the  sensors'  field  of 
view  is  called  the  Probability  of  Detection  (Pp)*  The  rate 
at  which  false  target  declarations  occur  is  termed  the  "False 
Alarm  Rate  (FAR)." 

One  may  correctly  Infer  from  the  statistical  nature  of 
the  two  principal  IRST  performance  measures  (FAR  and  Pp)  tha.*" 
the  IRST  detection  problem  is  Inherently  statistical. 

The  means  by  which  the  thre.shold  in  Fig.  1  is  estab¬ 
lished  will  not  bn  considered  until  Section  VI.  For  the 
purposes  of  the  present  discussion  ,  while  the  possibility  of 
time  variation  is  retained,  the  threshold  >’^(1)  assumed 
t  be  d e  t  e  rtii  I  n  i  .y  f;  i  c  .  ^  ^ 

''^^Consiatent  with  a  convention  of  random  process  theory, 
stochastic  quantities  are  assigned  capital  letters,  with 
sample  values  designated  by  the  oorrespi^nding  lower  case 
letters.  E.g..  a  sample  function  of  the  random  process  Y(t) 
is  de.signatod  y(t). 
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Analysis  oi  threshold  comparison  receivers  with 
stochastic  thresholds  Yjj(t)  is  also  of  interest,  particularly 
when  Y^(t)  is  derived  as  a  function  of  the  photocurrent  X(t). 
The  "adaptive  threshold"  radar ^ ^ and  sonar^ ^^^receivers 
are  of  this  variety.  An  IRST  adaptive  threshold  processor 
analysis  is  presented  in  Section  VI. 


1 1 .  OUTL INE  OF  RESU LTS 


1 ' 


A .  Objectives 

The  objective  of  the  following  section  is  to  outline  u  new 
method  for  predicting  the  expected  number  of  times  that  the 
photocurrent  from  an  infrared  sensor  will  cross  a  time-varying 
threshold,  as  the  sensor  scans  a  non-uniform  scene.  The  method 
is  based  on  the  development  in  Section  III.B  of  an  expression 
for  the  "threshold  crossing  rate"  ihj(t)  of  the  random  current 
Y(t)  (cf.  Fig.  1).  The  rate  ttj  is  Integrated  to  derive  the 
probable  number  of  threshold  crossings  mj(0,T)  during  a  time 
Interval  T^^; 


mj(0,T) 


/  ift-(t)dt 

T  ^ 
o 


Cl) 


where  is  the  time  interval  (t|  T/2. 

The  expected  number  of  threshold  crossings  raj(0,T)  on  the  time 
Interval  may  be  related  to  the  usual  False  Alarm 

Rate  (FAR)  and  Probability  of  Detection  (Pj^)  statistics  typically 
used  to  characterize  the  performance  of  IRST  devices,  as  dis- 
cus.sed  in  Appendix  C. 

An  expression  for  the  threshold  crossing  rate  tfij(t)  will  be 
presented  in  Section  II. C.  An  approximate  asymptotic  formulation 
for  mj  is  presented  in  Section  II. B.  The  approximate  e.xpression 
for  mj  has  sufficient  accuracy  in  many  cases  and  eliminates  the 
need  to  perform  the  integral  of  eq.  (1).  The  asymptotic  result 
given  here  Is  derived  in  Section  III.D. 

It  is  assumed  that  the  average  photocurrent  «^f  the 

current  X(t)  (cf.  Fig.  1)  is  a  known  function  of  time. 


6 


I 


As  discussed  in  Appendix  A,  is  a  function  of  the 

aperture  irradiance,  the  Modulation  Transfer  Function,  and 
several  other  sensor  parameters  of  the  system  to  be 
modeled  (the  "model  sensor"). 


7 


B.  Asymptotic  Formulation  for  mj 

The  basic  quantities  needed  to  perform  the  calculation 
for  the  expected  number  of  threshold  crossings  are 

a)  The  time-varying  average  value  m  (t)  of  the  random 

X 

current  X(t)  (cf.  Fig.  1),  and 

b)  The  transfer  function  H(f)  of  the  post-detector 

filtf'r  In  Fig.  1. 

Equivalently,  the  Fourier  transform  of  H(f)  can  be 
specified  in  place  of  H(f): 

h(t)-r{H(f)}  (2) 

The  function  h(t)  is  called  the  "impulse  response"  of  the 
post-detector  filter,  and  }  denotes  the  Fourier  transform 
of  the  bracketed  quantity. 

As  discussed  in  Appendix  A,  H(f)  is  presumed  to  include 
the  frequency  characteristic  of  the  photodetector.  Thus, 
the  photodetector  in  Fig.  1  is  "idealized,"  in  the  sense 
that  it  is  presumed  to  have  a  perfect  all-pass  electrical 
frequency  characteristic,  the  frequency-dependent  part  of 
the  detector  responslvity  having  been  lumped  together  with  H(f). 

As  shown  in  Appendix  A,  the  average  value  m  (t)  of  the 

X 

random  current  X(t)  in  Fig.  1  is  derivable  in  terms  of: 

a)  high  spatial  resolution  radiometric  imagery  of  the 
infrared  scene,  in  the  optical  passband  of  interest 

b)  the  postulated  modulation  transfer  function  (MTF) 
of  the  optical  train  of  the  model  system 

c)  the  quantum  efficiency  and  physical  dimensions  of 
the  focal  plane  detectors,  and 

d)  the  focal  plane  scan  velocity. 
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The  basic  equation  for  obtaining  is  eq  (A-14). 

The  asymptotic  expression  for  the  expected  number  of 
threshold  crossings  m^  can  be  expressed  In  terms  of  the 
quantities  m  a  ,  m, ,  and  a  ,  each  of  which  can  be  derived  as 

y .  y 

a  function  of  m„(t).  The  definitions  ofm  o  ,  and  c  , 

X  y  I  y  ^ 

are  as  follows: 


my(t)  -  E{Y(t)} 

(3a) 

m,(t)  =  E{Y(t)>  -  A 
z  y 

(3b) 

Oy^(t)  -  El[Y(t)  -  my(tif 

(4) 

-  E{[Y(t)  -  m^(t)]^} 

(5) 

Thus  ,my  is  the  mean  value  of  the  output  Y(t)  of  the  post 

detector  filter  (cf.  Fig.  1);  is  the  variance  of  Y(t), 

and  a  ^  is  the  variance  of  Y,*  The 
z 

expectation  operators  in  eqs .  (3)-(5)  denote  ensemble 
averages  over  the  photon  fluctuation  noise.  It  is  assumed 
throughout  the  analysts  that  photon  noise  is  the  dominant 
noise  mechanism;  in  this  case,  the  detector  is  said  to  be 
operating  in  the  ’’background  limited  performance"  (BLIP) 
regime . 


♦In  the  following  analysis  a  dotted  quantity  denotes  the 
time  derivative  of  the  corresponding  undotted  quantity. 
Also,  E{-}  denotes  the  expected  value  of  the  bracketed 
quantity.  9 


As  shown  in  Appendix  B,  m  o  ,  and  a  bear  the 

y  >  y  z 

following  functional  relationships  to  m  (t): 

iDyCt)  -  h(  t)  #  ^ 

Oy^(t)  -  e  h^(t)  ©  inj^(t)  (7) 

a/(t)  -  e  [h(t)]^  »  m^(t)  (8) 

where  h(t)  is  given  by  eq.  (2),  e  is  the  electronic  charge, 
and  #  denotes  the  process  of  convolution: 

QQ 

a(t)  ®  b(t)  ■  S  a(u  )  b(t  -y)dv  (9) 

^OP 

The  asymptotic  formulation  for  mj  requires  knowledge 
of  m  ,  0  ,  and  o  at  a  discrete  set  of  times.  In  order  to 
illustrate  the  two  types  of  time  points  that  have  special 
significance,  reference  is  made  to  Fig.  2.  Fig.  2a  is  an 
illustrative  example  showing  the  mean  value  m  (t)  of  the 

X 

random  current  X(t),  as  a  function  of  time.  The  corresponding 
ra„(t),  as  obtained  from  eq.  (6),  might  appear  as  in  Fig.  2b. 

As  shown  in  Fig.  2b,  the  threshold  function  y^is  taken  as 
a  constant  in  order  to  simplify  the  discussion  of  this 
section;  generalization  of 
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the  present  results  to  time-varying  thresholds  is  performed 
in  succeeding  sections. 

With  reference  to  Fig.  2b,  times  t^  and  tg  are  of 
special  significance  in  performing  the  calculation  for 
mj(0,T) . 

At  time  t„  the  mean  current  crosses  the  threshold 
with  a  positive  slope. 

As  shown  in  Section  III.D.,  the  crossing  counter 
fflj  is  incremented  by  unity  at  time  tg.  It  may  be  seen 
that  the  following  equations  arp  satisfied  at  time  tgi 


(12) 


rflyCtg)  >  0 

When  there  are  no  solutions  to  eq.  (12)  during  the 
time  interval  of  interest,  the  value  of  mj  is  generally 
dominated  by  contributions  to  the  integral  in  eq.  (1)  that 
occur  near  local  maxima  in  m  (t),  such  as  the  maximum  at 

y 

time  in  Fig.  2b.  The  asymptotic  expression  for  the  incre¬ 
ment  in  ra.  associated  with  the  local  maximum  of  at  t, 
u  y  1 

can  be  written  as; 

raj(tj^,T)  =  ‘^j(^l)  •  "Stj^  (13) 
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where  mj(tj^)  ta  the  threshold  crossing  rate  at  time  t^,  and 

6t,  is  the  effective  time  Interval  during  which  m  (t)  remains 

at  its  closest  proximity  to  the  threshold  y^.  The  time  t  In 

eq  (13)  is  relatively  unimportant;  the  choice  T>35t^  is 

generally  adequate  to  capture  the  principal  contribution 

to  fflj  that  occurs  in  the  neighborhood  of  t^. 

A  reasonable  approximation  to  ihj(t^)  valid  when  the 

time  fluctuations  in  m  (t)  and  ni„(t)  are  slow  compared  to  a 

X  y 

characteristic  fluctuation  time  for  h(t)  is  given  by: 


Aj(t^)  - 


(14) 


where  <?y(tj^)  and  obtained  from  eqs  (7)  and  (8), 

evaluated  at  time  t^,  and 


(16) 


Eq.  (14)  is  Rice's  original  equation  for  the  mean  threshold 

crossing  rate  of  a  stationary  Gaussian  noise. 

The  effective  proximity  time  6tj^  in  eq  (13)  is  given  by 

pTT  -j  J  (16) 

-  L’STETrmrT. 


where  'u(tp  is  the  second  derivative  of  u(t)  with  respect  to 
time,  evaluated  at  t  •  t^^. 


C .  The  Exact  Expression  for  riij 

Eq  (13),  while  simple  to  evaluate,  is  not  always  an 
accurate  approximation  for  mj.  As  borne  out  in  the  numerical 
examples  of  Section  V.  ,  eq.  (13)  becomes  progressively  less 
accurate  as; 

a)  the  target/background  contrast  decreases,  and 

b)  the  width  of  the  peak  at  (cf.  Fig.  2a)  decreases. 

The  implication  of  a)  and  b)  is  that  eq.  (13)  is  unsuitable 
for  analyzing  IRST  performance  in  the  regime  whore  the 
clutter  discrimination  problem  is  likely  to  be  most  serious. 
The  objective  of  this  section  is  to  present  the  complete 
expression,  eq.  (17),  for  the  mean  crossing  rate  thj 
derived  in  Section  III.C.  The  exact  expression  for  ihj 

is  not  subject  to  limitations  a)  and  b)  above.  In 
addition,  the  derivation  leading  to  eq.  (17)  is  used  to  guide 
the  derivation  in  Section  VI  of  a  rigorous  formulation 
for  adaptive  threshold  processors. 

The  exact  expression  for  the  mean  crossing  rate 
derived  in  Section  III.C.,  appropriate  for  non-stationary 
Gaussian  random  processes,  is 

,  2  * 
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Vbere  a„(t),  a,(t),  and  u(t)  are  given  by  equations  (7), 
y  z 

(8),  and  (15),  respectively.  The  functions  (|)(.)  and 
are  the  Gaussian  density  function  and  its  integral, 
respectively: 


2 

iji(x)  -  <2Tr)  exp  (-x^/2) 


X 

<t(x)  ■  / 


4>(z)  dz 

The  quantities  p  and  p^  in  eq  (17)  are  defined  by: 


(18) 

(19) 


P  -  [o-d-r  )  ] 


+  ru0^  -  u  = 


Po 

Pi 


(20) 


with  lUy  given  by  eq  (0),  and  u(.)  defined  as  the  unit 
step  function: 

(0  x<0 

y(x)  -  ji  x>0  (21) 


Eq  (17)  is  applicable  even  for  time-varying  thresholds  yQ(t), 
as  may  be  Inferred  by  noting  the  presence  of  :^Q(t)  in  eq.  (20). 

The  reason  for  including  y^Ct)  and  j^^Ct)  as  arguments 
of  rftj  in  eq.  (17)  is  explained  in  Section  II. E. 

I 

The  only  still-undefined  quantity  in  eqs.  (17)  and  (20)  is 
r(t),  defined  now  as: 


r(t)  s  CygjCt)  /  (Oy 


(22) 


where  C  is  the  cross  covariance  function  of  Y(t)  and  y(t) 
yz 

(cf.  eqs  (3)-(5)): 


Cy2  5  E{[y(t)  -  my(t)]  [Y(t) 
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m^(t)J  } 


(23) 


As  sbovm  In  Appendix  B  , 


Cy^(t)  -  OyCt)  Oy(t)  (24) 

Thus,  from  eqs  (22)  and  (24) 

r(t)  -  (25) 

It  may  be  noted  that  eq.  (17)  was  derived  subject  to 
the  assumption  that  the  photocurrent  y(t)  (of.  Fig.  1)  is 
a  non-stationary  Gaussian  process.  The  justification  for 
this  assumption  and  the  constraints  it  imposes  on  the 
validity  of  eq.  ri?)  are  discussed  in  Appendix  P. 

f'or  the  special  case  of  uniform  scenes  and  stationary 
thresholds,  it  follows  from  eqs.  (6)  and  (7)  that 

"“b  “  <>y  "  ^o  "  °  ^26) 

From  eqs.  (20),  (25),  and  (26) 

r  -  p  «  pj  -  0  (27) 

From  (17)  and  (27) 

-4 

which  is  identical  to  eq.  (14).  Thus,  eq .  (17)  reduces 
correctly  to  Rice's  original  equation  (14),  appropriate 
to  the  case  of  stationary  Gaussian  noise. 
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Eq.  (17)  appears  almost  Identical  to  an  equation 
'  ^  14  ^ 

derived  by  Cramer  and  Leadbettor.  '  The  relationship 
between  eq.  (17)  and  the  corresponding  equation  In 
Ref.[!l4l  is  discussed  In  Appendix  D, 
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D .  The  Character  of  the  Results 

The  character  of  the  results  obtained  with  eq ,  (17) 
is  illustrated  with  the  aid  of  Figs.  2-5.  Once  again,  the 
fact  that  the  process  Y(t)  (cf.  Fig.  1)  ha.s  a  time-varying 
mean  m  (t)  qualifies  Y(t)  generically  as  a  "non-stationary" 
process.  Also,  while  the  theory  developed  in  this  paper  is 
general  enough  to  accomodate  both  time-varying  and  stochastic 
threshold  functions,  the  present  example  is  simplified  by 
the  assumption  of  a  constant  threshold  y^  (as  depicted  in 
Pig.  2b). 

Figures  3  and  4  depict,  respectively,  the  crossing  rate 

thj  and  expected  number  of  threshold  crossings  mj  corresponding 

to  the  mean  current  m..(t)  and  constant  threshold  y  of 

y  o 

Fig.  2b  •  The  following  general  obuorvations  are  made  with 
respect  to  figures  2-4. 

1.  As  seen  in  Fig.  3,  the  function  thj(t)  is  appreciably 

non-zero  over  only  a  very  small  fraction  of  the  total 
time  interval.  The  pulse-like  waveforms  in  Fig.  3 

are  much  narrower  than  the  corresponding  peaks  in  Fig.  2. 

2.  If  solutions  exist  to  the  equation 

my(t.j)  -  y^  (29) 

subject  to  the  constraint  that 

ih(t)>0  ,  at  t  -  (30) 

than  the  character  of  my  is  determined  almost  entirely 

<J 

by  the  be  havior  of  thj  in  the  near-neighborhood  of  the 

times  t  . 
s 

*By  convention,  dotted  quantities  denote  time  derivatives 
of  the  corresponding  undotted  quantities. 
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For  each  solution  of  equation  (29)  that  satisfies 
constraint  (30),  i.e.,  for  each  positive  slope 
threshold  crossing  by  the  mean  current  m  (t),  the 

i7 

expected  number  of  crossings  tRj  Is  Incremented  by 
unity.  Time  in  figures  2-4  is  an  example  of 
this  effect. 

3.  If  m  (t)  lies  below  y  on  the  interval  (tJ<T/2, 

Jr  w 

i.e.,  if  no  solutions  exist  to  eq.  (29),  then  the 
character  of  mj  Is  determined  almost  entirely  by 
the  behavior  of  ihj  in  the  near-neighborhood  of  the 
local  maxima  of  .  Time  t^  in  Figures  2-4 
illustrates  this  effect. 

4.  The  precise  form  of  m  (t)  is  typically  unimportant. 

y 

For  example,  the  m  (t)  of  Fig.  5  might  yield  effectively 

y 

the  same  m,  as  the  m  (t)  of  Fig.  2,  subject  only  to  the 
0  y 

constraints  that  the  peak  value  and  curvature  of 
m(t)  at  time  and  the  values  of  tj^,  t2,  and  y^  are 
the  same  in  figures  2  and  5. 


It  is  assumed  throughout  much  of  the  analysis  that  the 
comparator  devico  in  Fig.  1  recognizes  only  those 
threshold  crossings  for  which  the  following  two 
inequalities  are  observed; 


i'(t)  >  ^^(t) 

(31) 

i'(t)  >  0 

(32) 

This  assumption  on  the  nature  of  the  comparator  device 
has  been  incorporated  into  the  analysis  in  such  a  way 
that  the  analytical  expression  for  rtij  is  insensitive 
to  downcrossings ,  such  as  that  at  time  t.^  in  Fig.  2. 

In  fact,  similar  formulations  tor  the  crossing  rate 
may  bo  developed  subject  to  different  assumptions  than 
(31)  and  (32),  o.g.,  that  the  comparator  device  is  sen¬ 
sitive  only  to  downcrosaings ,  or  both  upcrossings  and 
downcrossings .  In  general,  it  might  be  desirable  to 
incorporate  both  upcrossing  and  downcrossing  sensitive 
detectors  in  a  single  IRST  receiver,  in  order  to  develop 
an  e^3timato^  of  object  size. 

It  should  be  recognized  that  a  hardware  implementation 
of  the  crossing  detectors  analyzed  in  this  paper  must 
first  estimate  quantities  y(  t)  and  y(t)  before  tests 
(31)  and  (32)  cun  be  performed. 
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E,  Adaptive  Threshold  Processors 

The  results  described  thus  far  pertain  to  calculations 
of  the  crossing  rate  performance  of  processors  configured  as 
in  Fig.  1,  having  time-variable,  deterministic,  threshold 
functions  y^Ct).  However,  as  discussed  in  Section  VI  and 
Appendix  E,  potential  perfoi-mance  advantages  accrue  to  the 
use  of  threshold-establishing  circuitry  that  adapts  the  value 
of  the  threshold  to  the  prevailing  background  conditions. 

These  "adaptive  threshold"  processors,  such  as  that  shown  in 
Fig.  20,  give  rise  to  threshold  functions  that  are  not  only 
time-varying,  but  are  also  Inherently  stochastic, 

By  choosing  a  suitably  large  value  for  the  delay  time 
Tjj  in  Fig,  20,  the  filtered  current  process  Y(t)  and  the 
threshold  process  Y^(t)  are  rendered  statistically  independent. 
As  shown  in  Section  VI. B.,  it  follows  that: 

riijCt)  •  Ey  ^  (205) 

with  rti.(t|y  ,y  )  given  by  eq .  (17),  The  expectation  operator 

vl  O  O 

E  .  in  eq.  (205)  is  defined  by: 


where  f  ,  is  the  joint  density  function  of  the  threshold 
^0^0 

process  Y  and  it's  time  derivative  Y 
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Happily,  the  crossing  rate  expression  for  time-varying 
deterministic  thresholds,  eq .  (17),  is  found  to  be  also 
applicable  to  the  stochastic  threshold-crossing  rate,  eq.  (205). 

Eq.  (205)  is  further  developed  in  Section  VI. C,  and 
the  following  approximate  result  is  obtained: 


m 


f  ^  1 

^(t)  -  mj(t|  my^.m^^)  [l  +  \ay  / 


(225) 


2  2 

where  and  are  calculated  from  eqs.  (7)  and  (8), 

respectively.  Expressions  for  m  ,  a  ,  and  a,.  ,  are  obtained 

^0^0  ^o 

by  straightforward  analogy  to  eqs.  (6)-(8): 


(t)  -  K  h^(t-T^)  ^  m^(t)  (189) 

(t)  =  eK  hj  (190) 

0^  (t)  =  eK  j^Ro(t-T^)j  ^®mjj(t)  (191) 


The  gT,in  constant  K,  time  delay  T^,  and  impulse  response  h^(t), 
in  eqs.  (189)-(191)  are  design  parameters  of  the  adaptive 
threshold  scheme  shown  in  Fig.  20.  Expressions  for  the  terms 
and  Eg  in  eq.  (225)  appear  as  eqs.  (226)  and  (227)  in 
Section  VI . C . 
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As  discussed  in  Section  VI. D.,  the  threshold  variance 
2  2 

factors  a  and  o.  in  eq.  (225)  represent  a  "False  Alarm 
^o  ^o 

Penalty"  that  causes  the  false  alarm  performance  of 
adaptive  threshold  schemes  to  be  inherently  inferior  to 
fixed  threshold  schemes,  against  uniform  background  scenes 
It  Is  also  noted  that: 

a)  The  analysis  of  Section  VI  Is  applicable  to  a 
considerable  variety  of  adaptive  threshold  processors 
Including  the  structures  depicted  in  Figs.  22-24. 

b)  A  list  of  important  design  parameter  trade-off 
issues  is  identified  at  the  end  of  Section  VI, D.  The 
only  important  obstacle  to  performing  these  analyses 
is  the  lack  of  high  spatial  resolution,  radiometric, 
infrared  background  imagery. 
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Ill .  ANALYSIS 

A.  Introduction 

The  objective  of  section  III.B.  is  to  derive  an 
expression  for  the  expected  value  of  the  number  of  times 
a  nonstationary  noise  Y(t)  crosses  a  nonstationary  stochastic 
threshold  Y^Ct).  The  main  result  of  section  III.B.,  eq.  (43)i 
is  not  used  in  its  complete  generality  until  the  adaptive 
threshold  processer  analysis  of  Section  VI.  Eq.  (42)  .is 
particularized  at  the  end  of  Section  III.B.  to  the  case  of 
deterministic,  time-varying  thresholds  y^Ct). 

The  crossing  rate  equation  that  results,  eq,  (48), 
requires  knowledge  of  the  joint  density  function  fy^Cy,^) 
of  the  current  Y(t)  and  its  time  derivative  Y(t). 

As  discussed  in  Section  III. C..  the  current  Y(t)  is  a 

non-stationary , Gaussian  process.  It  follows  that  Y(t)  and 

« 

Y(t)  are  jointly  Gaussian  processes.  The  joint  density 
fy^  can  then  be  expressed  in  terras  of  five  characteristic 
statistics . 

The  explicit  expression  for  the  jointly  Gaussian  fy^ 
is  then  uf ed  with  the  general  crossing  rate  expression  (48) 
derived  in  Section  III.B.  to  derive  an  explicit  expression 
for  the  crossing  rate,  eq.  (70).  The  expressions  for  the 
five  important  current  statistics  given  in  Secti.on  III  .C.  are 
derived  in  Appendix  B.  The  complete  expression  for  the 
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expected  auraber  of  crossings,  eq.  (72),  generally  requires 
the  numerical  Integration  of  a  somewhat  complicated 
Integrand.  Section  III.D.  is  devoted  to  deriving  some 
approximations  to  eq.  (72)  that  don't  require  numerical 
integration,  and  that  have  greater  intuitive  appeal. 
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B.  A  Basic  Equation  for  Curve  Crossing  Rates 

The  first  objective  of  this  section  is  to  derive  eq. 

(39)  for  the  average  zero-crossing  rate  of  a  random  process 
G(t). 

The  process  G(t)  Is  not  required  to  be  either  stationary 
or  Gaussian.  Fig.  6  depicts  a  sample  function  g(t)  of  the 
process  G(t).  It  is  noted  that  g(t)  has  a  zero-crossing  at 
t  ime  t  q: 

BCt^)  -  0  (33) 

It  follows  from  (33)  and  the  well-known  properties  of  the 
Dirac  delta-function,  6(.),  that 


to+A 


1  ■  /  6(g)  dg  -  / 

-"1 


fi[8(t)]  lft(t)l  dt 


(34) 


where  ,  eg,  and  A  are  defined  in  Fig.  6. 

It  follows  from  (34)  that  the  total  number  of  zero  crossings 
of  the  sample  function  g(t)  on  the  time  interval  lt|  <  T/2 
is  given  by : 


J  “  /  'S(g)  lg(t)Idt 

^o 


(35) 


In  writing  (35)  it  is  assumed  that  g(t)  does  not  have  a  zero¬ 
crossing  either  at  t  ■  i  T/2.  Defining  nij(0,T)  as  the 
expected  value  of  J  on  the  time  interval  : 
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raj(0,T) 


E{J} 


(36) 


and  operating  on  oq.  (35)  with  the  expectation  operator 
E{ . } ,  one  obtains 

MjCO.T)  -  /  dt  E(<S(8r)|ji(t)  I }  s  /  ihj(t)dt  (37) 

'^o 

with  the  interval  defined  by  |t|  <  T/2. 

Prom  (37) 

Ihj(t)  -  E{6(g:)|4|}  -  /di[  /dg  (g.i)  {5(g)  14|  }  (38) 

where  (gi4)  is  the  Joint  probability  density  function  of 
the  process  G  and  its  time  derivative  (j. 

The  first  basic  result  of  this  section  is  the  following 
expression  for  the  zero-crossing  rate  thj  : 

Aj(t)  »  ;l6|  (0,fe)  d4  (39) 

Eq.  (39)  is  well-known;  however,  most  references  to 
(39)  appear  to  impose  a  statlonarity  requirement  on  G  that 
is  not  actually  necessary.  The  applicability  of  (39)  to 
nonstationary  processes  appears  to  have  first  been  recognized 
by  Cramer  and  Leadbetter^ 

The  domain  of  integration  in  eq.  (39)  is  a  matter  of 
some  interest.  If  one  wishes  to  calculate  only  the  expected 
number  of  positive  slope  zero-crossings,  i.e. ,  the  expected 
number  of  times  that  both 

G(t^)  -  0  (40) 

and 

G(t^)  >  0 
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(41) 


are  both  satisfied  on  the  Interval  |t|  <  T/2i,  the  lower  and 
upper  limits  of  integration  In  (39)  should  be  chosen  as  0  and 
<»,  respectively.  The  resulting  expression  for  mj(0,T): 

00 

mj(0,T)  -  /  dt  {/  d4  \&\tQ^(0  ,it))  (42) 

o  ^ 

does  not  include  zero-crossings  of  the  type  depicted  in 
Fig.  6,  for  which  i  <  0.  Apparently,  thj  is  sensitive  only 
to  the  "right  type"  of  zero  crossing,  as  defined  by  the 
limits  of  integration  in  eq.  (39). 

The  process  G(t)  is  now  assumed  to  be  formed  as  the 
difference  of  two  stochastic  processes  Y(t)  and  Y^jCt)  : 

G(t)  =  Y(t)  -  Y^(t)  (43) 

where  Y^(t)  is  referred  to  as  the  "threshold  process". 

Without  making  any  assumptions  with  respect  to  the  statistics 
of  Y(t)  and  y^Ct)  (e.g.,  each  process  may  be  both  non¬ 
stationary  and  non-Gaussian)  it  follows  from  (43)  that: 


;/d;dn  f 


+  n,c 


n) 


(44) 


The  proof  of  (44)  is  straight-forward  (cf.  ref.  12113  ,  p.l31). 

Eqs.  (42)  and  (44)  are  the  basis  of  the  adaptive  threshold 
processer  analysis  of  Section  VI  .  For  the  present,  it 
suffices  to  specialize  (44)  to  the  case  of  non-random 
thresholds,  for  which  the  fourth-order  Joint  density  ^ 

becomes 

(.z,k  +  n.c.n)  -  fy^Cc.d:  +  n)  s(y^-z)  6(S'Q-n)  (44a) 


From  (44a)  and  (44) 

-  'y^<yo*'k  +  ^0> 

Prom  (39)  and  (48) 

ilij(t)  -  /|4|  f  y^(y^,i  +  5'o)dj:  (46) 

Changing  variables; 

z  s  4  +  (47) 

(46)may  be  written  as 

iftj(t)  -  /|z-^o|fyj^(yQ.z)dz  (48) 


The  quantity  rfij  ,  previously  interpreted  as  the  zero-crossing 
rate  of  the  process  G  (cf.  eq.  (39)),  is  now  interpreted  as 
the  threshold  crossing  rate  of  the  process  Y(t). 
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The  dotnalQ  of  integration  of  eq.  (48)  is  deflnr«d  by 
assuming  that  the  hardware  implementation  of  the  threshold 
counter  is  sensitive  only  to  those  crossings  for  which  the 
following  two  inequalities  are  satisfied  simultaneously: 


o 

A 

(49) 

V 

O 

(50) 

For  the  example  of  Fig.  7  this  particular  type  of  detector 
would  increment  its  "crossing  counter"  at  time  t^,  but,  by 
design,  would  fall  to  Indicate  threshold  crossings  at  times 
tg,  tg,  and  t^. 

It  follows  from  eqs.  (48)  -  (50)  that 

ifaj(t)  -  !(<»)  -  (51) 

Where  u(.)  Is  given  by  oq.  (21),  and 

X 

I(x)  S  /  (a-:^^)  fy^(yjj,z)dz  (52) 

The  crossing  rate  thj  could  Just  as  well  have  been 
formulated  subject  to  different  constraints  than  (49)  and 
(60).  For  example,  the  hardware  Implementation  of  the 
threshold  counter  may  be  sensitive  to  all  crossings  for  which 

:>(t)  >  0 

regardless  of  the  relative  magnitudes  of  J’(t)  and  ^^(t). 
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In  this  case  eq.  (48)  Is  still  subject  to  the  constraint 
of  eq.  (50),  but  the  constraint  of  eq.  (49)  is  no  longer 
in  force.  It  may  then  be  shown  that  the  following 
expression  obtains,  instead  of  eq.  (51); 

'»j(t)  -  I(-)  -  (53) 

with  I(.)  given  by  (52). 

When  <  0i  satisfaction  of  (50)  assures  satisfaction 
of  (49).  One  then  expects  and  ibj  to  be  equal.  Inspec¬ 
tion  of  eqs.  (51)  and  (53)  shows  that  this  is  indeed  the  case: 

ifaj(t)  -  rfij(t)  -  I(«)  ,  ^^(t)  <  0 

However,  when  >  0  condition  (49)  la  more  restrictive 
than  (60).  One  then  expects  ihj  to  be  greater  than  rtij. 

From  (51)  and  (53): 

rfij(t)  -  ttj(t)  -  I(S-q)  ,  :^Q(t)  >  0 

Since  1(9^)  is  inherently  negative  (cf.  (52)),  it  follows 
that  >  thj,  as  expected,  when  >  0- 

With  respect  to  Fig.  7,  it  can  be  shown  that  the  integral 
iTij  of  eq.  (63)  is  incremented  by  unity  both  at  time  t^^ 
and  tg.  By  comparison,  mj  is  incremented  by  unity  only  at 
time  tj^i  ttj  does  not  "count"  the  crossing  time  tg,  for 
which  condition  (49)  is  violated.  Neither  mj  nor  Wj  register.s 
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a  crossing  at  times  tg  and  in  Fig.  7,  since  condition 
(50)  is  violated  at  times  tg  and  t^. 

Plots  oi  m(t)  and  ffl(t)  corresponding  to  Fig.  7  are 
depicted  in  Figs.  8a  and  8b,  respectively. 

The  analysis  that  follows  in  Sections  III  -  V  Is 
based  on  development  of  eq.  (51)  rather  than  eq,  (53), 
l.e,,  it  Is  assumed  that  the  signal  processer  hardware 
being  modelled  requires  the  satisfaction  of  both  (49)  and 
(50)  before  it  will  register  a  threshold  crossing. 
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c. 


N0N-3TATI0NARY  GAUSSIAN  PROCESSES 


Further  development  of  eq.  (48)  for  the  crossing 
rate  ibj  Is  contingent  on  obtaining  a  suitable  expression 
for  the  Joint  density  f^^.  The  objectives  of  this  section 

are: 

1)  to  evaluate  eq.  (48)  for  the  particular  case  of 
a  bivariate  Gaussian  density,  eq.  (C5),  and 

2)  to  ihtroduoe  the  expressions  derived  in  Appendix  B 
for  the  five  basic  statistics  that  appear  as 
parameters  in  eq.  (56). 

The  Justification  for  assuming  a  Gaussian  distribution 
for  f  (and  hence  for  f.  and  f  .  as  well)  is  given  in 

y  y  y  y 

Appendix  F. 

The  approximation; 

Fj,(y)  “  «  (  Oy  ) 

I 

( 22 ) 

has  been  considered  by  Gilbert  and  Pollack,  'where  '!>(.) 
is  the  Gaussian  distribution  function,  eq.  (05),  ra  and 
a-  are  defined  by  eqs.  (56)  and  (58),  and  'F  (y)  is  the 

J'  i/ 

distribution  function  of  the  process  y(t): 

F  (y)  -  J-y  f„(y')dy' 

^  ..OQ  ^ 
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A  summary  of  Gilbert  and  Pollacks'  results  is  given  in 
Ref.  (15). 

The  problem  addressed  in  Appendix  F  is  different, 
and  much  simpler,  than  the  problem  of  Ref.  (22).  It  is 
shown  in  Appendix  F  that  the  relative  error  in  thj  is 
approximately  equal  to  the  relative  error  c  in  the  density 
function  of  Y( t ) ; 

fyCy)  “  ♦  ("5^)  o  +  ") 

where  (|i(.)  is  the  Gaussian  density  function,  eq,  (04), 

(13) 

The  Edgeworth  series  expansion^  ■^of  f  (y)  provides  a 

¥ 

Simple  and  easily  evaluated  expression  for  the  relative 

error  £.  Sample  calculations  described  in  Appendix  F 

show  that  £  is  negligibly  small  for  typical  system  and 

background  parameter  values. 

Based  on  the  analysis  of  Appendix  F,  the  joint  density 

f  .  is  now  a.ssumed  to  have  the  following  bivariate  normal 

y  y 

form : 


oi  ru^  +  -  2ruv'| 

fy^Cy^j.z)  -  {2TTayajj(l-r  )}  exp  {-  L  2(l-r^)  ■*  ^ 

where : 


niy  =  E{Y} 

(B6) 

■  E{Y) 

(57) 

dj  S  El(Y-nly)^) 

(58) 

5  E{(Y-m^)^) 

(59) 

=  E{(Y-niy)(Y-m^))  2  Cy.^ 

(60) 

u  2  (y„-my)/dy 

(61) 

V  = 

(62) 

34 


It  may  be  shown  from  eqs  (52)  and  (55)  that 


I(x) 


(1-r^)  <j>(u){  GKPq)  -  +  P;i^C'f(P^-I>(Px)^  > 


(63) 


where  a  ,  a  ,  r,  and  u,  are  given  by  eqs  (58),  (59),  (60), 

y  z 

and  (61),  respectively.  The  functions  (()(.)  and  4'(.)  in 
(63)  are  defined  as  follows: 


(Ji(x)  =  (2ti)“^  exp  (-x‘'/2)  (64) 

X 

^(z)dz  (65) 

The  quantities  p^,  Pj^ ,  and  Pj^  in  (63)  are  defined  as: 

p^  =  f(l-r^)*  (m^j  +  ruo^)  (66) 

p^  =  pl-r^)^  (m^  +  (67) 

Px  “  C(^”^  )  (m  +  rua  -x)  (68) 

Z  Z 


From  eq.  (68),  (64),  and  (65) 

lim  <ti(p  )  -l.im  (I)(p  )  »  0 


(69) 
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From  eqa.  (51),  (63),  and  (69): 

o  i 

ttj(t)  -  Voy/(l-r“)  (J)(u)  {•|i(p)  +  Pjl'(p)} 
wh'ire  by  definition 


(70) 


(71) 


with  p^  and  given  by  aqa.  (66)  and  (67). 

The  expected  number  of  threshold  crossings  on  the  time 
interval  •  defined  by  [tj  <  T/2,  is  given  by: 

?n,,(0,T)  -  E  {J}  «  /  ih.(t)dt  (72) 

V  U 

O 

with  rhj(t)  given  by  eq.  (70).  In  order  to  complete  the  formal 

calculationul  procedure  for  mj,  it  is  necessary  to  specify 

values  for  the  five  functions  m„(t),  m„(t),  a„(t),  a_(t), 

y  z  y  z 

and  r(t',  defined  by  eqs.  (56)  -  (60).  It  is  shown  in 
Appendix  B  that  for  background-limited  (BLIP)  detection, 


my(t)  ■  h(t)  ^ 

o^ct)  -  e  h^(t)  W  m,^(t)  (74) 

o^(t)  -  e  Cri(t)J^  «  m^{t)  (75) 

m^(t)  “  rhy(t)  (73a) 

r(t)  «  ^y(t)/a,^(t)  (76) 


whoi'o  e  la  the  oloutronio  ohiirge,  h(t)  is  the  Impulse 
roaponso  of  the  post-dotoctor  riltiu*  (of  Fig.  1), 
is  the  moan  value  of  the  ourront  X( t)  (of.  FIr.  1),  and 
♦li  is  the  convolution  operator: 

f(t)  •'♦K(t)  -_l  f(t-x)g(.x)dx  (77) 

The  relationship  of  ihj^(t)  to  the  radlunoo  of  the 
scene  under  observation  and  the  optical  parameters  of 
the  IRST  sensor  is  disouasod  in  Appendix  A,  and  expressed 
quantitatively  by  eq.  (A-l-l).  Eqs.  (A-14),  (70),  and 
(72)-(7G),  taken  together,  represent  a  proscription  for 
calculating  the  performance  of  a  time-variable  threshold 
IRST  sensor  against  a  non-uniform  .scene, 

In  principal,  it  remains  only  to  model,  measure,  or 
otherwise  specify  the  radiance  distributions  ^^f  typical 
targets  and  background  scones.  Thl.s  aspect  of  the  IRST 
modeling  problem  Is  not  considered  In  this  paper. 

It  may  bo  noted  that  eq.  (70)  appears  almost  Identical 
to  an  equation  .in  Ref.  (1*1).  The  r«'lation.shtp  holwt'en 
eq .  (70)  and  t.ht*  cori’e.spondliig  equathm  in  Ucl’,  (M)  1  .s 
discussed  In  Appendix  D. 


D.  Asymptotic  Approximations 

1 .  Introduction 

Evaluation  of  eq .  (72)  for  the  mean  number  of  threshold 
crossings  TOj  generally  requires  a  numerical  integration  of 
tftj,  as  given  by  eq.  (70).  It  is  now  shown  that  a  relatively 
simple  approximation  for  mj  can  be  developed  that  eliminates 
the  need  to  perform  this  numerical  integration,  and  has 
greater  intuitive  appeal  than  the  original  formulation. 

Some  numerical  examples  comparing  the  numerical  integral  of 
eq.  (70)  with  the  approximate  result  developed  here  are 
presented  in  part  V  of  this  paper.  Eq.  (72),  with  eq.  (70) 
substituted  for  thj,  is  now  written  as; 


mT(0,T)  -  /  a(t)e~'^^^^dt 

o 

where  is  the  time  interval  |t)  <T/2,  and  where 


(78) 


b(t)  5  u^(t)/2 


a( 


t) .  (S 


)  (4>(p)  +  P;j^  'l’(P)}  (2TT) 


-i 


(79) 


(80) 


Experience  with  the  numerical  evaluation  of  eq.  (78) 
has  led  to  the  conclusion  that  ihj  is  appreciably  non-zero 
only  over  very  brief  intervals  of  time  (of.,  for  example, 

Fig.  3).  Since  the  principal  contribution  to  the  integral 
mj  of  liij  accrues  over  those  same  short  intervals  (cf.  Fig.  d), 
oq.  (78)  may  be  written  as: 
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(81) 


mj(O.T)  «  j:  mj 
n 


where 


wan  m 

n 


(82) 


Where  T„  Is  the  time  interval  1 t-t„ I  <  t„/2. 

The  interval  is  presumed  to  be  large  enough  for  the 
integral  in  eq.  (82)  to  capture  the  principal  contribution 
to  raj,(0,T)  made  by  the  n^^  term  in  the  summation  of  eq.  (81). 

Experience  with  the  numerical  evaluation  of  eq .  (78) 
has  also  shown  that  the  rapid  time  variation  of  ihj  in  the 
neighborhood  of  the  discrete  times  t^^  is  mainly  due  to  the 
exponential  factor,  the  function  tt(t)  varying  much  more  slowly 
with  time.  Eq .  (82)  is  thus  approximated  as: 


m,(t„,r„)  a  a„I_ 
u  n  n  n  n 


(83) 


where  a^^  =  ‘^(In).  and 


I  =  /  e'^^^^dt 


( 84 ) 


where  T„  is  the  time  Interval  |t-t  I  <  t„/! 
n  '  n '  n 


The  usual  first-order  saddle  point  method  is  now  used 
to  obtain  approximations  for  1^^.  The  saddle  points  t^  aro 
found  by  solving  for  the  roots  of  the  equation: 
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(85) 


B(tn)  -  u(tjj)  Ci(t^  )  -  0 


tq.  (85)  has  two  distinctly  different  types  of  solutions, 
uorrosponding  to  solutions  of  the  following  two  equations; 


Mt^)  -  0 


(86) 


u.;t„)  -  0 


(87) 


Each  root  of  eq.  (85)  gives  rise  to  a  discrete  contribu¬ 
tion  to  the  approximate  formulation  for  mj,  the  form  of  which 
depends  on  whether  the  root  is  a  solution  to  eq .  (86)  or 


eq.  (87). 


Solutions  of  eq .  (86),  such  as  time  tj^  in  Fig.  2b, 
will  be  referred  to  as  "saddle  points."  Solutions  of 
eq.  (87),  such  as  time  tg  in  Fig.  2b,  will  be  referred  to 
as  "mean  crossing  times." 

The  next  objective  is  the  development  of  approximations 
to  the  integral  1^^  of  eq .  (84). 

2 .  Saddle  Points 

The  function  u(t)  in  eq .  (79)  is  now  given  a 
Taylor  expansion  about  a  time  t^^; 


u(t)  % 


a,^(t-tn)“/2 


(88) 


Assuming  that  t^  is  a  "saddle  point",  i.e.,  a  solution  of 
eq .  (86),  eq.  (88)  becomes: 


u(  t ) 


U  ( t  - 
n' 


•5 

)“/2 


(89) 
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From  eqs.  (79)  and  (89),  dropping  the  fourth-order  term: 


b(t)  -  u^/2  +  u^u^(t-t^)2/2 


(90) 


From  eqs.  (84),  (90),  and  (79): 

e  /2 

'v,  “  ..2 

I„  “  exp(-b„)  f  0xp(-u^  u  t  /2)dt 
n  «  -e^/2  “ 


(91) 


where  =  b(tj^). 

Uaklng  use  of  the  fact  that 

«3cp  (-irx^)dx  =•  1 

and  assuming  that  In  (91)  Is  chosen  sufficiently  large 
eq.  (91)  can  be  approximated  as 

-  (27r/u^  u^)^  exp  (-b^) 

From  (83)  and  (93) 

Comparison  of  eqs.  (78)  and  (72)  shows  that 
ifaj(t)  =  a(t)  exp  |^-b(t 
With  the  definition: 

i  (2,./u„  li„)‘ 


(92) 


(93) 


(94) 


(95) 


(96) 
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It  follows  from  (94)-(96)  that: 


From  eqs.  (81)  and  (97) 


"‘j'O.’')  -  I  “j^b)  «‘b 


It  Is  relatively  straightforward  to  show  that: 

a(t)  5  d/dt  {(yo-my)/Oy}  --o^(l-r^)*  p^/o 


Thus,  from  (86)  and  (99) 


Pl(t^)  -  0 
From  (80)  and  (100) 

“  (o^/Oy)(l-r^)*  (t(p^)  (2tt)"* 


From  (100),  (86),  (67),  and  (71) 


=0  .  <  0 


Pi  -  0 


at  t  ■  t^. 

n 


From  (95),  (101),  and  (79) 


2oi 


ihj(tjj)  -  {(2tt)"*  (Oz/^7  )  ^(Ujj)}  {[2ff(l-r^)0®  ^(Pjj)} 


(97) 


(98) 


(89) 


(100) 


(101) 


(102) 


(103) 
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The  desired  approximation  to  mj  is  given  by  eq. 

(98),  with  and  rh.(t„)  given  by  eqs.  (96)  and 
n  j  n 

( 103  )  .respectively. 

The  first  bracketed  quantity  on  the  left-hand-side 
of  eq .  (103)  is  identical  to  Rice's  equation  for  the 
threshold  crossing  rate  of  stationary  Gaussian  noise.  The 
second  bracketed  quantity  in  eq.  (103)  is  always  less  than 
unity.  For  constant  threshold  processors  (^-^  ■  0)  the 
second  bracketed  quantity  is  equal  to  ( 1-r  )'*,  which  is 
nearly  unity  whenever  the  fluctuations  in  the  mean  current 
m  (t)  are  slow  compared  to  the  impulse  response  h(t)  of 
the  post-detector  filter. 

For  constant  threshold  processors,  the  quantity  in 
eqs.  (96)  and  (98)  can  be  interpreted  as  the  effective  time 
interval  during  which  the  filtered  mean  current  m  (t)  remains 
in  the  neighborhood  of  its  peak. 

It  should  be  noted  that  the  approximation  (98)  may 
be  incomplete,  since  the  discrete  terms  in  the  summation 
all  correspond  to  roots  of  eq.  (86).  As  discussed  in  the 
following  section,  roots  of  eq.  (87),  assuming  such  roots 
exist,  each  make  an  additional  discrete  contribution  to  mj. 

3 .  Mean-Crossing  Times 

Instead  of  eq .  (86),  it  is  now  assumed  that  eq .  (85)  is 
satisfied  because 

U(t„)  -  -  0 


A.3 


(104) 


The  saddle  points  obtained  as  solutions  of  eq.  (104)  correspond 
to  times  when  the  mean  current  «>y(t)  crosses  the  threshold 
y^j(t).  Thus,  in  the  vicinity  of  t^, 
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Subject  to  assumption  (109),  it  follows  from  eq.  (108) 
that 


♦(p„) 


0 

m 


and 


jKrhy) 

iKifiy-^o) 


:^o  < 


where  U( . )  denotes  the  unit  step  function 


0 

li<x)  -  1 


X  <  0 
X  >  0 


Eq.  (Ill)  can  also  be  written  as 

't(Pn)  -  M<thy)  U(rhy  -  p^) 

Noting  that 

Vi(*y  “  Pq)  Bgn  (rtiy  -  y^)  =  u(rtiy  -  P^) 

it  follows  from  eqs.  (107),  (110),  (113),  and  (114),  that 

subject  to  assumption  (109),  with  given  by  (108). 

It  follows  from  eq.  (115)  that  If  the  mean  current 
i«y(t)  crosses  the  threshold  y^^(t)  at  time  t^^,  then  the 
crossing  count  itj  will  be  incremented  by  unity  if  the 
following  two  conditions  are  both  satisfied: 


(110) 


(111) 


(112) 


(113) 


(114) 


(115) 


45 


1)  th«  Llope  of  TOyCt)  at  time  (l.e.,  ih(t^))  must 
be  positive,  and 

2)  the  slope  of  m  (t)  must  be  greater  than  that  of 
y^^(t)  at  time  t^^ 

—assuming  that  the  slopes  of  ra  (t)  and  y_Ct)  at 

JF  W 

time  aren't  nearly  the  same,  i.e.,  that  ®y(t)  and  yQ(t) 
don't  run  nearly  parallel  in  the  neighborhood  of  their 
crossing  point  t^.  (This  last  assumption  is  expressed 
quantitatively  by  inequality  (109).) 

Thus,  the  asymptotic  result  for  mean-crossing  times 
reflects  the  limitations  on  the  domain  of  integration  in 
eq,  (48)  Imposed  by  inequalities  (49)  and  (50). 
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IV.  INFINITE  CONTRAST  POINT  TARGETS 


A.  Introduction 

The  mean  cro.<3sing  rate  formalism  developed  in 
Section  III.C  requires  adherence  to  the  following  pre¬ 
scription; 

1.  The  mean  current  m  (t)  must  be  specified.  (It 
Is  shown  in  Appendix  A  how  ni^(t)  may  be  calculated 
In  terms  of  the  radiance  distribution  of  the  scene.) 

2.  The  post-detector  filter  impulse  response  h(t) 
must  be  specified. 

3.  Eqa .  (73)-(76)  must  be  evaluated  for  the  five 

functions  my(t),  ra^jCt),  and  r(t). 

4.  The  mean  crossing  rate  thj(t)  is  obtained  from 
eq.  (70),  and  integrated  over  time  to  obtain  the 
expected  number  of  crossings. 

The  difficult  calculational  problem  generally  posed 
by  step  3  above  can  be  circumvented  by  making  the  following 
choice  for  m  (t ) : 

m^it)  -'Xq  6(t-t^)  (116) 

Eq.  (116)  is  the  limiting  case  of  an  Infinite  contrast 
point  target,  assuming  an  ideal  optical  MTF,  and  point 
detectors  in  the  focal  plane. 
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\ 


II 


The  final  reauxt  of  this  section 
for  in^(0,T)  that  is  formally  identical 

U 

sion  'or  the  probability  of  detection 


will  be  an  expression 
to''the  usual  expres- 

p  (16). 

D 


P 


(117) 


Where  <f(.)  is  the  Gaussian  distribution  function,  and 
tp  is  the  time  a*-,  which  the  post-detector  filter  output 
current  assumes  its  peak  value. 

iThe  application  of  eq.  (117  to  calculations  was 
(16) 

previously  derived  subject  to  the  assumption  of  internal 
amplifier-noise 'limited  operation,  for  which  the  filtered 
current  Y(t)  is  a  stationary  random  process.  The  surprising 
result  is  now  obtained  that  eq .  (117)  is  also  applicable  to 
Pp  calculations  for  BLIP  sensors,  for  which  the  current  Y(t) 
is  a  nonstationary  random  process.  The  only  important  dis¬ 
tinction  between  the  calculational  procedures  for  for 
amplifier-noise-limited  or  background-limited  operations  lies 
in  the  means  for  obtaining  o  (t  ).  For  BLIP  sensors,  ci  (t  ) 

y  r  J  Ir 

is  calculated  from  eq .  (146). 
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B. 


Analysis 


Eqs.  (116)  and  (73)-(76)lead  immediately  to  the  following 
equations  for  the  five  basic  current  statistics: 


my(t)  » 

x^hCt-to) 

(118) 

a/(t)  - 

(119) 

Xq  e  [Mt-t^) 

2 

(  120) 

ihy(t)  =» 

Xq  h(t-t^) 

(121) 

r(  0  cTyCt)  “ 

h(t-t^) 

(122) 
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i 

i 

Prom  oqs.  (119)  and  (120): 

qy(t)  «  !h(t-t^^)|  -  (x^o)^  h(t-t^)  sgn  [^(t-t^^)]  (123) 

a2(t)  -  lh(t“t^)i  -  (x^^e)^  h(t-t^)  sKn  [^(t-t^)]  (124) 

where  the  function  agn  (.)  denotes  the  sign  of  Its  argument. 

Prom  eqs.  (122)-(124): 

•  r(t)  -  agn  [h(t-t^P 

Thus,  from  (126): 

ir(t)i  -  1  (126) 

Unity  correlation  coefficient,  eq.  (126),  invalidates 
the  bivariate  density  eq,  (56),  (Note  the  factor  ( 1-r  )“  in 
eq.  (55).)  Thus,  a  new  expression  for  f  muat  be  found  that 

y  « 

is  valid  for  the  ci'se  r  ■  1,  and  used  to  evaluate  eq .  (52). 

It  may  bo  shown  that  eq.  (126)  implies  the  following 

\ 

relationship  betwoeji  Y(t)  3  Z(t)  and  Y(t): 

Z(t)  -  a(t)  Y(t)  (127) 

where 

a(t)  2  {h(t-t^^)/h(t-t^j)}  (128) 
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It  follows,  in  turn,  that: 


.z)  ■  i^iy)  -  ay) 

Where  i5(.)  is  the  Dirac  delta  function,  a  is  given  by 
(128),  and 

with 


<j)(x)  =  (2Tr)“^  exp  (■-x^/2) 


and 


u(t)  = 

From  eqs .  (52)  and  (129) 


y^(t)  -  niy(t)j 


I(x)  -  (ay^  -  Sr^)  S,^(y^)  I m( ay^ )-M (ay^-x) } 

with  the  unit  step  function  u(.)  defined  by  eq .  (112). 
Assuming  a  constant  threshold  processor, 


-  0 


it  follows  from  eqs.  (51)  and  (133)  that 


rtij  (t)  =■  !('»)  *=  ay^  f^  (y^,)  P(ay^^) 


(129) 


(130  ) 


(131) 


(132) 


(133) 


(131) 


(135) 
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From  eqs.  (72),  (135),  and  (130)-(132), 


mj(0,T)  -  /  dt  ay^  (t.(u) 

With  the  following  change  of  variable 

I  • 

<1  =  ho®)^ 


(136) 


(137) 


eq.  (136)  may  be  written  as 
g(T/2)  , 


ij(0,T)  -  ^  U  [B(t-tQ)yQ/h(t-t^)j  sgn  [h(t-tQ)j  (|)(q)dq 


(138) 


It  Is  now  assumed  that  the  threshold  Is  positive 


y  >0 
•^o 


(139) 


and  that  h(t)  has  the  form  shown  In  Fig.  9..  Eq.  (138)  then 
may  be  written  as; 


q(tj^) 

ra,(0,T)  -  -  / 

q(t„) 


({)(q)dq 


(140) 


From  eq.  (137) 


Yo-^,  h(0) 


«  lim  l(x^e)*  h(0) 
h(0)-*-0'^ 


(141) 
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V.  SAMPLE  CALCULATIONS 


A.  A  Model  Filter  and  Model  Background 

In  ordor  to  calculate  the  threshold  crossing  rate  rfij,  it 

is  generally  necessary  to  first  numerically  evaluate 

2  2 

eqs.  (73)-(76)  for  the  functions  m  (t),  a  (t),  o  (t), 

Jr  jr  w 

m  (t),  andr(t),  since  analytical  evaluation  of  eqs.  (73)  - 
(76)  is  generally  not  possible. 

In  order  to  obviate  the  need  for  numerical  evaluation 
of  (73)  -  (76),  it  is  assumed  in  this  section  that  the 
impulse  response  h(t}  and  mean  current  m  (t)  have  the 
following  forma : 

h(t)  ■  (S/t^)  exp  {-fT(t/t^)^}  cos  Snf^t 

mj^(t)  «  Xq  {1  exp  [-ir(t/T)^j} 

Eq.  (149)  is  a  two-parameter  family  of  functions,  with 
parameters  and  f^.  Eq.  (150)  is  a  three-parameter  family, 
with  parameters  x^,  x^,  and  t.  Eqs.  (149)  and  (150)  are 
particularly  convenient,  as  they  allow  analytical  evaluation 
of  eqs.  (73)  -  (76),  and  have  a  large  enough  number  of 
parameters  to  illustrate  many  of  the  effects  that  one  would 
expect  to  observe  with  more  general  h(t)  and  m  (t). 

<ra 

Taking  the  Iburier  transfomti  of  eq.  (149)  results  in  the 
following  expression  for  the  transfer  function  H(f): 


(149) 

(150) 
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(151) 


H(f)  -  F{h(t)}  -  {exp  +  exp  [-TT(f+f^\2j) 

Representative  plots  of  eqs .  (150)  and  (151)  are  given  in 
Figs.  (11)  and  (12). 

It  is  seen  from  Fig.  11  that  the  parameters  x^,  x^^,  and 
t  of  ra  (t)  correspond  respectively  to  a  constant  background 
brightness  level,  the  object/background  peak  contrast,  and 
the  size  of  the  non-uniformity.  It  follows  from  Fig. 12 
that  the  model  filter  has  a  bandpass  characteristic  centered 
at  the  frequency  f^.  From  eq.  (149),  the  parameter  t^  is 
a  damping  constant  for  the  filter.  Thus,  while  not  completely 
general,  eqs.  (149)  and  (150)  permit  the  exploration  of  a 
considerable  variety  of  possible  background/filter 
interaction  effects. 

Assuming  that, 

exp  (-ATTf^^t^^)  >>  1  (152) 

it  follows  from  eq.  (151)  that 

max  H(f)  ^  1  (153) 

f 

The  factor  (2/t^)  on  the  right-hand-side  of  eq.  (149)  was 
chosen  to  normalize  H(f)  to  a  peak  value  of  unity. 

The  only  major  deficiency  in  eq.  (151)  as  a  model  filter 
characteristic  is  the  inability  to  control  the  filters' 
low  frequency  roll-off  independently  of  the  center-frequency 
parameter  f^  and  noise  bandwidth  parameter  t^ . 

The  derivations  of  the  analytical  e.xpressions  that  result 
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from  evaluating  eqa.  (73)-(76)  with  the  assumed  functional 

forms  for  h(t)  and  m  (t)  are  available  from  the  author. 

Figs.  i3-17  are  illustrative  evaluations  of  m  (t), 

¥ 

2  2 

o„(t),  C„„(t),  and  a„(t),  assuming  the  following 
z  y  yz  z 

choices  for  the  five  model  background/filter  parameters; 

X  -  134.  namps  (154) 

o  ■ 


f  -  6.91  kHz  (157) 

o 

t  -  225  nsec.  (158) 

o 

The  rationale  behind  the  choice  of  parameter  values  (154)  - 
(158)  is  discussed  in  Appendix  G. 
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B.  Crossing  Count  Calculations 

In  this  section  the  value  obtained  for  tnj(0,T)  by 
integrating  eq.  (70)  numerically  will  be  compared  with 
the  value  obtained  for  mj(0,T)  by  three  approximate  methods. 
In  every  case,  it  will  be  assumed  that  h(t)  and  m^^Ct)  have 
the  forms  specified  by  eqs.  (149)  and  (150),  and  that  the 
threshold  is  independent  of  time. 

The  first  approximation  to  mj(0,T)  is  given  by: 

m^gCO.T)  -  /  jtiQg(t)dt  (159) 

where 


ihj^g(t)  s  <27t)i  ^160) 

and  is  the  interval  |tj  <  T/2. 

Comparing  eq<  (160)  and  eq.  (14),  it  is  seen  that  eq.  (159) 
is  simply  the  integral  of  Rice's  original  expression  for  the 
crossing  rates  of  stationary  processes,  transmuted  into  a 
function  of  time  through  the  use  of  time-variable  expres¬ 
sions  for  the  basic  current  statistics,  m  ,  a  and  a  . 

JT  JF  ^ 

It  is  generally  found  that  the  estimate 

mqs  “  "'J 
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Improves  monotoalcally  as  the  parameter  x  in  eq.  (150)  is 
made  larger.  Eq.  (159)  defines  what  will  be  referred  to  as 
the  "quasi-stationary"  (Q-S)  approximation  to  the  mean 
crossing  count  mj. 

The  second  approximation  to  mj(0,T)  that  will  be 
evaluated  in  this  section  is  the  asymptotic  approximation 
derived  in  Section  III.D: 

m^(0,T)  H  rt>j(tg)6tg  (162) 

The  saddle  point  t  is  obtained  front  eqs.  (100)  and  (67)  as 

s 

the  solution  of 

^  (163) 

since  it  is  assumed  that  ^'^(t)  «•  0  in  this  section  in  order 
to  simplify  the  numerical  examples.  It  has  been  found  that 
(163)  is  effectively  equivalent  to 

ra2,(tg)  •  0  (164) 

due  to  the  very  small  value  of  r(t„)  for  the  examples  that 
have  been  worked.  It  follows  from  eq.  (164)  and  inspection 
of  Fig. 14  that  t  ■  0  for  the  example  of  eqs.  (154)  - 

o 

(158).  The  same  saddle  point  condition  has  also  been  found 
to  hold  for  the  other  numerical  cases  discussed  in  this  section 
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t  -  0 

s 


(165) 


corresponding  to  the  time  of  closest  approach  of  niy(t)  to 
y^,  for  the  assumed  expressions  (149)  and  (150)  for  h(t)  and 
m  (t).  It  follows  from  the  above  discussion  and  the  analysis 
of  sectionlll.D  that  the  appropriate  asymptotic  approximation 
to  mj(OiT)  is  given  by 

ra^(O.T)  -  AqgCO)  6tg  (166) 

where  is  given  by  eq.  (160),  and  6tg  is  given  by 
eq.  (96) : 

fitg  -  {2Tr/u(0)  u(0)l^  (167) 

where 

u(0)  -  Oy(0  ) 

and 

11(0  )  =  d^u(t)/dt^ 

t=0 

The  third  and  final  approximation  to  inj(0,  T)  is  effectively 
a  modified  version  ofOenoud's  method^ ^ for  the  calculation 
of  probability  of  detection; 

nipg(0,T)  -  1  -  4'{u(0)}  (170) 


(168) 


(109) 
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where  ')(.)  is  the  Gaussian  distribution  function,  as  given 
by  eq.  (65). 

As  will  presently  bo  shown,  eq .  (170)  becomes 
increasingly  accurate  as  an  approximation  to  raj.; 

mpg(0,T)  «  mj(0,T)  (171) 

as  the  object  contrast  is  made  progressively  larger 
and  the  object  extent  t  is  made  progressively  smaller. 

This  is  not  at  all  surprising,  in  light  of  the  fact  that 
the  mean  threshold  count  for  infinite  contrast  point 
targets,  eq.  (144),  is  identical  to  eq.  (170).  The 
quantity  mpg  will  bo  called  the  "point  source  approximation" 
to  mj. 

Fig. 18  is  a  plot  of  the  base  ton  logarithms  of 

mj(0,Tp),  mpg(0,Tp),  and  m(j^(0,Tp),  as  a  function  of  the 

normalized  peak  threshold  level  u^ : 

P 

Up  2  u(0)  (172) 

where  u(0 )  Is  given  by  eq.  (168),  and  u.ssuming  the 
background  radiance  and  filter  parameter;.;  of  eqs.  (154)- 
(158). 
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As  discussed  In  Appendix  G,  parameters  (154) -(  158) 
correspond  to  the  following:  more  basic  parameter  choices; 


a  ■  1  mrad 

0 

(object  angular  subten.se) 

01  ^  ■  1  mrad 

(instantaneous  fiold-of- 

-view) 

Tp  -  1  .sec 

(frame  time) 

(173) 

H(0)  -  10“'^ 

(transfer  function  aero 

ordinate ) 

n  *  0 .  G 

(quantum  efficiency) 

(174) 

\  ■  4 . 0pm 

(optical  wavelength) 

o 

A.,„  ■  500  cm“ 
ap 

(aperture  area) 

-2 

Ljj  ■  0,14  mw  cm 

sr"^  (background  radiance) 

( 17.5) 

Specifying  ^  qualifies  the  object  described  by 

eq.  (150)  as  ''targf?t-lil<e"  in  character  (cf.  Appendix  G). 

It  is  seen  from  Fig.  18  that  the  point  source  appro.xi- 

mation  mpg  provides  an  e.xcellent  fit  to  the  crossing  count 

function  nij,  for  the  target-like  object  described  by 

eqs.  (154)-(15S).  However,  the  asymptotic  approximation 

m.  underestimates  m .  by  nearly  an  order  of  magnitude  for 
A  ».» 

the  target-like  object. 


61 


It  is  now  assumed  that  the  object  is  clutter-like, 
rather  than  target-like: 

»  5  0.  *  5  mrad. 

0  T 

instead  of  the  previously  assumed  values  *  1  mrad. 

It  follows  that 


T  ■  795  usee. 

replaces  eq.  (156),  The  other  background  and  filter 
parameter  values,  eqs,  (154),  (155),  (157),  and  (158)  are 
left  unchanged.  Fig.  19  is  a  plot  of  the  base  ten 
logarithms  of  the  corresponding  values  of  mj(0,Tp), 
mpg(0,Tp),  and  m^(0,Tp)  for  the  clutter-like  object. 

It  is  seen  from  Pig.  20  that  the  asymptotic  approxi¬ 
mation  is  effectively  equal  to  nij,  providing  a  much 
better  approximation  than  nipg  for  the  clutter-like  object. 

It  has  generally  been  found  for  a  fairly  wide  range  of 
background  and  filter  parameters  that 

a)  rupg  provides  an  excellent  approximation  to  m^ , 

and  is  superior  to  m^,  when  the  time  .scale  of  fluctua¬ 
tion  of  m  ( t )  is  comparable  to  that  of  h(t). 

b)  m^  provides  an  excellent  approximation  to  in^,  and 
is  superior  to  Wpg,  when  the  time  scale  of  fluctuation 
of  m  (t)  i.s  slower  than  that  of  h(t). 
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VI. 


Adaptive  Throshold  Pr(>o.Hasora 


•  Inti’oductioti 

The  basic  IRST  processor  under  consideration  up 

to  this  point  is  described  in  Section  I.B.  and  doplct«id 

in  Fig.  1,  It  has  been  shown  in  Section  V  that  the 

threshold  crossing  performance  of  the  receiver  in  Fig. 

1  is  strongly  dependent  on  the  threshold  function 

y^(t).  Thus,  it  is  highly  desirable  that  the  IRST 

processor  suppress  clutter-induced  threshold 

crossings  by  increasing  v  (t)  when  m  ( t )  is  "i^lutter- 

o 

like'' . 

Rather  than  allow  .v^^(t)  to  take  on  an  a  priori 
constant  or  functional  value,  it  is  necessary  to 
establish  the  thre.shold  by  some  means  that  "adapts" 
y^^(t)  to  the  prevailing  backgrouml  conditions. 

A  similar  type  of  signal  processing  problem  has 
been  addressed  in  the  radar^  ^and  .st)nar^ ^ 

Literatures.  A  candidate  adaptive  threshold  sclumu' 

(  '-Xj  \ 

adapted  from  the  earlier  work  Is  depictod  In 
Fig.  20. 


o4 
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The  block  with  transfer  function  exp(-j  2iTfT^) 
introduces  a  delay  of  seconds.  The  triangular-shaped 
block  in  Fig.  2Q  denotes  an  ideal  all-pass  amplifier  of 
gain  K.  It  is  seen  that  the  delay  time  T^,  gain  K,  and 
transfer  functions  H^(f)  and  H(f)  are  all  design  variables. 
Strategies  for  choosing  the  design  variables  in  order  to 
satisfy  particular  performance  requirements  or  optimization 
criteria  will  not  be  discussed.  The  remainder  of  this 
section  is  devoted  to  developing  a  formulation  for  the 
expected  number  of  threshold  crossings  for  IRST  receivers 
structured  as  in  Fig.  20  under  the  assumption  that  the 
design  variables  have  all  been  specified. 

Eq.  (225)  for  the  crossing  rate  rt)j  is  the  principal 
result  of  this  section. 

B.  Analysis 

The  starting  point  for  the  adaptive  threshold  analysis 
is  eq .  (39) : 

mj(t)  =  /|gl  fQ^(0.g)dfe  (176) 

where,  from  eq.  (43); 


(c  c,n) 


(177) 


Assuming  that  the  processes,  y ,  y^,  and  are  jointly 
Gaussian  (cf.  Appendix  F  for  justification),  their 
joint  density  can  be  expressed  in  terms<of  their  covariance 
matrix^  ^  A  .  The  matrix  ^  has  four  rows  and  four 
columns,  for  a  total  of  sixteen  elements.  Written  in 
partitioned  form: 


where  the  superscript  T  in  (178)  denotes  the  matrix 
transpose  operation.  The  sub-matrices  C,  and  in 
eq.  (178)  are  defined  as  follows: 


C  « 


ay(t) 


C^yCt.t) 


Cy^<t.t)) 


So  - 


0*  (t) 
^o 


'^O'^O 


^0^0 


(t,t) 


CT?  (t) 


Si 


^yoy(t,t) 


Cy„^(t,t)) 


c^yt,t) 


C^^^(t,t)( 


(179 


(180 


(181 


The  scalar  covariances  that  comprise  the  elements  of 
(179)-(181)  are  defined  by 


C^jj(ti,t2)  =  E{  A(tj^)-n^(tj^)  j^B(tj.)-mg(t2 ) 


'}  ( 182 


where  A  and  B  take  on  the  values  Y.  Y,  and  Y^,  as 


appropriate.  Also. 


where  h^(t)  is  the  Fourier  inverse  of  Taking  the 

expected  value  of  both  sides  of  (188): 

tn  (t)  -  K{h^(t-T.)  W  m^(t)}  (189) 

yo  o  a  X 

analagous  to  eq.  (D-30).  It  may  also  be  shown,  analagous 
to  (185)-(187)  that: 
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6y^^(t)  -  eK{mj^(t)  ®  h^Ct-Td)} 

(190) 

2 

(5^^(t)  -  eK{mj^(t)  »  [h^(t-Td)]^} 

(191) 

^  (t.t)  -Oy  (t)  (t) 

(192) 

It  remains  only  to  formulate  similar  expressions  for  the 
elements  of  in  order  to  complete  the  specification  of 
the  joint  density  ^  . 

The  analysis  this  requires  is  so  similar  to  Appendix  B 
that  it  will  only  be  outlined  here, 

It  may  be  shown  that: 

oa 

^yy^^l’'^2^  "  ^(tg-u)  (193) 

Also, 


^^1 

(194) 

V  <*r‘2>  - 

(195) 

(196) 
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From  (193)-(196); 
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From  eqs.  (176),  (177),  and  (203) 

m,(t)  -  IfdKdn  f  ^  (^  ,n){/d4lj!:|f  .  (<  ,4+n) }  (204) 

o  yy 

With  the  change  of  variable  z  ■■  i  n  •  eq.(204)  may  be 
written 


“j't’  *  <2055 

where,  by  definition, 

*j(tlyQ,^Q)  2  /  |z-^^lfy^(y^,z)dz  (206) 


Ey  ^  {.}  2  //dcdn  f„  *  (r.,n){.}  (207) 

o  ^o'o 

and  le  obtained  by  analogy  to  eq.  (66). 

Happily,  eq.  (206)  is  seen  to  be  identical  to  eq.  (48). 
Thus,  the  entire  analysis  of  sections  IV. B  and  IV. C, 
originally  applicable  only  for  deterministic  thresholds,  is 
now  found  to  be  directly  applicable  for  stochastic  thresholds 
as  well.  From  eqs.  (70)  and  (71).* 


ihj(t 


yo'^c^ 


(1-r^)^  <l-(u){4)(p)  +  Pj 


♦  (P)} 


(208) 


where 

jPo  3^0^° 

p  “  ['‘’y  ^^o^]  “ri 
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Eq.  (205)  reprssqnts  a  formal  means  of  oaloulatlng  the 
mean  threoholn  crossing  performance  of  the  adaptive 
threshold  processor  depicted  in  Fig.  20. 

However,  the  evaluation  of  eq .  (205)  appears  to 
present  some  significant  calculational  difficulties. 
These  difficmtlas  are  obviated  by  means  of  the  approxi¬ 
mate  method  of  evaluation  pursued  in  the  following 
section . 

C.  Approximations  for  ift. 


FdUomlng  ref.  (26),  p.  Ul,  *^(t|  la  given  a 
truncated  Taylor  expansion  about  the  values  y  -  m  and 
fo  ■  ° 


^  c/2 

where 


(210) 


A  2  *  (t 


=  rf)j(t|  .  ) 


(211) 


y,,  y  J 

o 


(212) 


:2i3) 
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The  following  notation  has  been  employed  in  eqs.  (212)  and 
(213): 


thj(tl.)  -  V^lyc^o^l 


y  ■  ni 
^o  y^. 

❖  -  rt> 

o  y. 


^0  "  ®y, 
^o  “  *y 


From  eqs.  (205)  and  (210) 


A.Ct)  «  iL  ^  (A)  +  E„  ^  {B}  +  J  E„  ^  {C} 
^  ^0*^0  ^0^0  ^0^0 


From  eqs.  (207)  and  (211)^(213) 


Ey  ,  U)  - 


O'  o 


EL  ^  {B}  -  0 


Ey  ^  {C} 


"’yo  ^yQyo  ^ 


o  ■'0‘'O 


From  (216)-(219): 

2  2  2  2 

rhj(t)  =;  (l+4(a  3  „  +  a .  3.  ^  )}  rfi  (t|  .) 

“  ■'O  ■’^O  O  '^O  ^0^0  “ 


(214) 

(215) 

(216) 

(217) 

(218) 

(219) 

(220) 
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Further  development  of  eq.  (220)  la  simplified  by 
employing  eq.  (D-6)  for  rhj,  rather  than  the  crossing  rate 
function  of  eq.  (208).  From  eqa .  (D-4)  and  (209): 

U(Pl)  +  't(P;]^)l 


(221) 


where 


Pi  “  (*y  - 


(222) 


Consistent  with  the  numerical  results  obtained  thus  far,  terms 
of  first  and  higher  order  in  r  have  been  dropped  in  writing 
eqs.  (221)  and  (222). 

From  eqs.  (220)  and  (221): 


-(i^)  [l  *  4  Mu); 


(223) 


y  "tn 
o  y^ 


where 


i  -  ‘KP^)  +  Pi4>(Pj^) 


(224) 


It  may  be  shown  from  (223)  and  (224)  that: 


rti.(t)  “  rtip,  (t  m 


.Ky  )  [i  (^)  +  (Sf ) 


(225) 
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where 


i  -  1) 

(226) 

•>'0 

and 

Cg  =  <j)(p^)/2i  I  (227) 

y  "m 

°  yo 

The  quantities  m  (t),  o  (t),  and  o.  (t)  in  eq.  (225)  are 
calculated  by  means  of  eqs,  (189),  (190),  and  (191),  respectively. 
Finally,  the  expected  number  of  threshold  crossings  in  a  time 
interval  may  be  obtained  as; 

mj.(0,T)  -  /  itij(t)dt  (228) 

o 

with  ibj  given  by  eq.  (225). 

Eq.  (225)  provides  the  basis  for  analyzing  a  much  broader 
range  of  possible  adaptive  threshold  schemes  than  Fig.  20  might 
suggest.  For  example,  straightforward  generalizations  of  eq.  (225) 
may  be  applied  to  the  structures  of  Figs.  22-24,  Further  numerical 
analysis  along  the  lines  of  Sec.V  for  various  model  backgrounds 
and  candidate  processor  structures  is  indicated,  but  has  not 
been  pursued  thus  far. 
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Some  qualitative  discussion  of  eq.  (225)  follows  in 
the  next  (and  concluding)  section  of  this  paper. 

D.  Discussion 

The  advantages  and  disadvantages  that  accrue  to  the  use 
of  an  Adaptive  Threshold  (AT)  processor  may  be  evaluated  by 
comparing  eq.  (225)  with  the  corresponding  crossing  rate  function 
for  Fixed  Threshold  (FT)  processors; 

rfaj,^(t)  -  rfi^L(tiy^,0)  (229) 


The  most  important  consideration  in  computing  the  False  Alarm 
Rate  (FAR)  of  FT  processors  is  the  factor: 


(t>(u) 


exp 


-4 


■)  ®  (230) 

that  enters  oq.  (229)  through  eq.  (221).  As  shown  in  the 

numerical  examples  of  Section  V,  constant  threshold  processors 

operating  against  non-uniform  scenes  are  most  susceptible  to 

false  alarms  in  the  immediate  neighborhood  of  the  local  maxima 

of  m  (t),  i.e.,  local  minima  of  u^(t).  It  follows  from  eq.  (230) 

that  the  false  alarm  rate  is  an  actremoly  sensitive  function 

of  m  ,  increasing  at  a  tremendous  rate  as  the  mean  current 
y 

niyCt)  approaches  y^. 

The  potential  advantage  of  AT  schemes  may  be  seen  by 
comparing  eq.  (230)  with  the  corresponding  factor  in  eq .  (225): 


(t)(u) 


(2tt)~^  exp 


V  “ra 

^o  y 


[- 


(m  -m  )^/2o  ^ 

y  y 


(231) 
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If  the  filter  H^Cf)  in  Fig.  20  can  be  chosen  such  that  m  (t) 

"tracks"  the  background- induced  variations  in  m  (t), 

2 

the  performance-destroying  local  minima  in  u  (t)  can  be 
eliminated. 

It  should  be  noted  that  the  adaptive  threshold  performance 
advantage  Just  described  is  only  realized  when  the  background 
scene  is  non-uniform.  The  performance  of  AT  processors  is 
generally  inferior  to  the  performance  of  FT  processors  when 
the  background  scene  is  uniform  and  of  known  brightness. 

In  this  case; 

FAR^^  -  FARp^Cl  t  FAP)  >  FAR^^  (232) 

where 

FAP  =  ^2 

The  threshold  variance  terms,  eq,  (233),  contribute  a  False 
Alarm  Penalty  (FAP)  whenever  an  AT  processor  is  used  when  it 
truly  isn't  needed. 

According  to  eq.  (232),  the  adaptive  threshold  false 
alarm  rate  (FAR^,p)  is  greater  (l.e.,  worse)  than  the  fixed 
threshold  false  alarm  rate  (FARp,j,),  assuming  that  the  adaptive 
threshold  gain  K  in  Fig. 20  has  been  adjusted  to  achieve 
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For  uniform  backgrounds,  it  Is  not  hard  to  show  that: 

(a  /a  -  (Af  /Af)  (230) 

und 

-  (Af^/Af)^  (237) 

where  Af  is  the  noise  bandwidth  of  H(f),  and  Af^  is  the  noise 
bandwidth  of  H^(f)  (cf.  Fig. 20  ).  It  follows  from  eqs.  (233), 
(236),  and  (237),  that: 

FAP  ^  (Af^/Af)  {e^  +  (Af^/Af)^  e (238) 

Generally,  Af^^  should  be  chosen  smaller  than  Af  to  minimize 
the  false  alarm  penalty,  eq .  (238),  and  to  prevent  a  too- 
rupid  threshold  response  that  would  tend  to  suppress  target- 
induced  threshold  crossings.  On  the  other  hand,  Af^  should 
be  chosen  large  enough  to  allow  the  threshold  to  accurately 
follow  most  of  the  structure  in  the  background  scene.  Clearly, 
the  choice  for  Af^  involves  degrading  system  performance 
against  uniform  backgrounds  for  the  sake  of  improved  performance 
against  non-uniform  backgrounds. 


Eq.  (237)  is  derived  by  assuming  a  rectangular-shaped  H(f), 
having  an  upper  cut-off  frequency  ,  and  a  noise  bandwidth 
Af  f^. 
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It  appears  likely  that  a  more  favorable  trade-off 
cQuld  be  achieved  with  the  receiver  structure  shown  in 
Fig.  22,  both  from  the  standpoint  of 

a)  decreasing  the  false  alarm  penalty,  eq .  (233),  and 

b)  improving  the  background  tracking  properties  of  m  (t). 
Eq.  (225)  and  the  entire  analysis  of  the  preceding  section  is 
easily  adapted  to  the  structures  of  Figs.  22-24.  The  false 
alarm  penalty,  eq.  (233),  decreases  roughly  as  (2N)”^  for  the 
detector  of  Fig.  23  ,  The  improvement  in  uniform  background 
performance  thus  obtained  for  large  values  of  N  is  gained  at 

the  expense  of  degraded  performance  against  cluttered  scenes, 
as  compared  to  detectors  with  small  values  of  N.  The  good 
background-tracking  capability  of  the  structure  in  Fig.  22 
combined  with  the  low  false  alarm  penalty  of  the  structure 
of  Fig.  23  can  be  obtained  by  employing  a  two-dimensional 
detector  array  with  'Hme-Cblay  and  Integration  (TDI)  logic. 

In  order  to  put  this  discussion  on  a  concrete  quantitative 
basis,  particular  background  and  target  radiance  distributions 
must  be  chosen,  and  the  mean  current  n'jj.(t)  calculated  by  means 
of  eq.  (A-14).  The  target  detection  and  clutter  rejection 
capabilities  of  a  given  candidate  adaptive  threshold  processor 
can  than  be  analyzed  by  means  of  eq.  (225).  Intercomparisons  of 
the  numerical  results  thus  obtained  for  a  variety  of  different 
processor  structures  should  than  allow  quantitative  conclusions 
to  be  drawn  concerning  such  issues  as; 
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a)  the  performance  penalty  caused  by  failing  to  match 
the  sensor's  Instantaneous  f ield-of-view  to  the  angular 
size  of  the  target . 

b)  the  potential  performance  advantages  of  time-delay 
and  integration  (TDI) . 

c)  the  best  value  of  N,  and  the  desireability  of  having 
different  transfer  functions  H^(f)  for  each  of  the  2N 
taps  in  the  tapped  delay  line  structure  of  Fig.  24. 

d)  the  advantages  that  may  be  gained  by  employing  simple 
two-dimensional  threshold  processing,  in  which  the  "target 
signal"  Y(t)  and  threshold  function  Y^Ct)  are  derived 
from  detectors  scanning  at  different  elevations. 

e)  the  possible  advantages  of  including  power-law  devices, 
or  other  instantaneous  non-linear  devices  in  the 
threshold-establishing  circuitry.  It  is  noted  that  the 
propagation  of  mean  values  and  variances  through  non¬ 
linear  processing  elements  may  be  treated  by  relatively 

straightforward  means. Thus,  the  calculation  of  m  , 

2  " 
and  0  for  use  in  eq.  (225)  is  not  difficult,  for  a 

specified  non-linear  relationship  between  Y^(t)  and  X(t). 

The  only  important  obstacle  to  performing  analyses  of  the 

kind  described  above  is  the  lack  of  high  spatial  resolution, 

radiometric,  infrared  background  imagery. 
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Fipure  I 


A  basic  threshold  comparison  receiver.  The  photodetector 

in  this  figure  is  "idealized",  in  the  sense  that  it  is  presumed 

to  have  a  perfect  all-pass  electrical  frequency  characteristic; 

the  frequency-dependent  part  of  the  detector  responsivity  is 

lumped  together  with  the  transfer  function  of  th.e  post- 

detector  filter  to  obtain  H(f).  A  "target  declaration"  is 

made  whenever  the  filtered  current  Y(t)  e.xceeds  the  thre.shcUd 

level  Y  (t). 
o  ' 
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Part  a)  of  this  figure  is  an  Illustrative  example  showing 
the  mean  value  m  (t)  of  the  random  current  X(t),  as  a  function 
of  time.  A  .small  "target-like"  (i.e.,  short-duration)  dis¬ 
turbance  occurring  at  time  t^  is  superposed  on  a  larger  ampli¬ 
tude,  more  slowly  varying  background.  The’ corresponding 
mean  current  m^Ct)  at  the  output  of  the  post-detector  filter 
(cf.  Fig.  1)  is  shown  in  part  b)  of  this  figure.  The  slowly- 
varying  background  part  of  is  greatly  attenuated.  By 

contrast,  the  target  pulse  is  attenuated  very  little,  since 
its  bandwidth  is  matched  to  that  of  the  post-detector  filter. 

A  constant  threshold  of  amplitude  y^  appears  as  a  hori¬ 
zontal  line  in  Fig.  2b),  The  time  t^,  where  the  mean  current 
m^Ct)  crosses  the  threshold  with  a  positive  slope,  is  referred 


to 


as  a  "mean-crossing  time"  in  the  text.  The  time  t, .  where 


7o  - 


my(t)|  has  a  local  minimum,  is  called  a  "saddle  point" 


in  the  text . 
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Figure  4 

”.'he  expected  number  of  threshold  crossings  nij  is  shown 

as  a  function  of  time,  obtained  as  the  integral  of  the 

crossing  rate  function  rtij  depicted  in  Fig.  3.  The  small 

Jump  in  mj  at  t^  is  associated  with  the  local  maximum  in 

ra  (t)  at  t,  shown  in  Fig.  2.  The  unity  increase  in  m^  at 

y  i  d 

time  tg  also  has  its  corresponding  origin  in  Fig.  2.  The 
absence  of  a  similar  effect  at  t^  (cf.  Fig.  2b)  is  due  to 
the  definition  of  rh.  as  a  positive-slope  cro.ssing  rate, 
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Figure  5 


The  average  current  tn  ( t )  depicted  in  this  figure  yields 

J 

nearly  the  same  curve  for  the  expected  number  of  threshold 

crossings  m.(0,t)  (cf.  Fig.  4)  as  the  in.  tt)  in  Fig.  2.  It  is 

u  y 

assumed  only  that  the  peak  value  and  curvature  of  trty(t)  u-t  time 
tj^  and  the  values  of  t^,  and  are  the  same  in  Figs.  2  and  5 
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Figure  6 

A  sample  function  g(t)  of  the  random  process  G(t)  is 
depicted  us  a  function  of  time.  The  particular  sample  function 


chosen  has  a  zero-crossing  at  time  t  •  The  time  A  is  chosen 

o 

small  enough  such  that  the  time  interval  |t  -  t^J  A  brackets 
no  zeros  of  gCt)  other  than  that  at  t  . 


y^(t)  and  n>yCt)  are  an  illustrative  threshold  function  and 
mean  current,  respectively.  When  the  domain  of  Integration 
in  eq.  (48)  is  chosen  to  satisfy  constraints  (49)  and 
(50),  the  expected  number  of  threshold  crossing.s  iTij  is 
incremented  at  time  tj^ ,  but  not  at  times  t2i  ,  or  t^^. 
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Figure  8 


Figures  8a  and  8ta  depict  the  time  development  of  the 
expected  number  of  threshold  crossings  calculated  by  means 
of  equations  (51)  and  (53),  respectively,  for  the  threshold 


y^(t)  and  mean  current  "'y('t)  shown  in  Fig.  7. 
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Ftjjuro  11 

An  1 1  Luatrat ivi?  plot  of  tho  model  baokffround  ourront 
la  shown  as  a  function  of  time.  The  parametric  form  of  the 
model  if>,^(t)  i.s  p:iven  by  oq .  (150).  The  piiramt?ters  ,  .x^,  and 
V,  correspond  respectively  t<5  a  oon.stunt  back|j:round  level,  the 
obJect/backKround  peak  contrast,  and  the  aiae  of  an  "object." 
in  the  scene.  DepeudinK  on  whether  r  i.s  larKO  or  small  compared 
to  the  .system  dwell  time,  the  model  current  •n,j(t)  la  representa¬ 
tive  of  a  "c  Lutter-L  ike"  or  a  "tar^^et-like"  object,  re.spectively 


Fl^ur*'  12 

A  ri'prtjaontut  I  VO  phit  ol’  the  modol  Tlltor  rransfor 
function  H(f1  Is  shown  as  u  function  of  the  ('loctrUral 
freipioncy.  The  puriuitotrlc  form  of  Hlf)  Is  irlvon  by  »'q . 
(151).  Tlie  purumotor  f^^  is  the  contor  froquoncy  of  tho 
bnndpusvS  chursictoi-lstio  .  Tho  purnmotor  l.s  a  docay  t  lino 
for  the  flUairs'  ImpaUso  rosponso  tof.  oq .  MP).  and  Is 
directly  related  to  tho  noise  bund\vldtl\  of  tiio  filter  by 
oq ,  (G-Tj).  The  cero-ord  Inatt*  H(0)  of  Hlf)  is  dotormlnod 
by  the  product  of  f^^  and  t^^  (of.  oq,  0— I ) ,  and  is  q-oin'rally 
a  very  small  number,  Tho  function  11(f)  Is  normallciv'd  to  a 
peak  value  of  unity. 
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An  lllustnil.  Ivo  ovaluatlon  of  ploftuHi  tis  u 

function  of  time,  for*  th»'  model  biio kK^'ound  ami  filter 
puriuneters  of  cqH.  ( 154  ) -( 158  )  ,  The  assumed  parametric 
fortn.s  for  the  model  filter  impulse  respon.st»  IH  t )  and 
the  model  hackKi'ound  pllven  by  eqs.  t  Ml))  aiui 

(150).  respectively,  The  mean  current  m^^lt)  plotted  in 
this  fixture  Is  defined  by  «'q .  (.50),  and  evaluated  in  terms 
of  h(  t )  and  ni,.(t)  by  means  of  eq .  (73). 
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A  plot  of  m.(t)  2  tti  (t)  is  shown  as  a  function  of  tlnu', 
y  y 

for  the  model  buckBround  and  filter  parameters  of  eqs.  (154)- 
(158).  This  curve  is  the  time  derivative  of  the  curve  in 
Figure  13. 


TIME  (msec) 

Figure  15 

9 

A  plot  oi  Oy“(t)  is  shown  as  a  function  of  time,  for 

the  model  parameters  of  eqs.  (154)-(158).  The  variance 
o 

Cy  (t)  of  the  filtered  current  Y(t)  is  defined  by  eqs.  (5 
and  evaluated  in  terms  of  h(t)  and  m. (t)  bv  means  of  eq. 
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-0.2  L. 

Figure  li) 

A  plot  ot  CynCt)  "  C^,^,(t)  Is  shown  as  a  function  of  time, 
for  the  model  parumetors  of  oqs .  il54)-(l58).  The  quantity 
Cy,^  is  the  orosH-covar  iance  of  the  rilt»'red  ourrent  Y(t)  and 
its  time  derivative  (cf.  oq ,  00),  and  is  k^valuated  hy  means  of 
oq,  (H-23). 
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FiKuro  17 

O 

A  plot  of  is  shown  for  t  lu>  iiiodo  L  punimotors  of 

*4 

•1 

otps.  (  1 3-1 )  -  ( 158 )  .  Tho  variance  o  “(.t)  of  i  ho  time  derivative 
of  the  filtered  current  is  defined  by  eq .  (59).  and  evaluated 
by  meaius  of  oq .  (7S). 
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NONMALiaO  THMtSHOLD 

Flguro  18 

The  base  ten  logarithm  of  the  expected  number  of  threshold 
crosalngu  mj  is  plotted  us  a  function  of  the  normaiiiied 
threshold  level,  for  the  filtered  current  statistics  of  Figs, 

13-17.  The  solid  curve  Is  the  quantity  mj .  obtained  by  inte¬ 
grating  eq.  (70)  numerically.  The  dotted  curve  is  the 
asymptotic  approximation  to  nij  given  by  eq .  (166).  The  dashed 
curve  tiipg  i.s  the  point  source  approximation  to  calculated 
from  eq .  (170).  The  two  approximations  to  mj  offer  calcula- 
tional  advantages,  in  that  neither  m^  nor  mpg  require  numerical 
integration  for  their  calculation.  The  basic  model  baclcground  and 
filter  parameters  used  in  calculating  these  curves  are  given  by 
eqs.  (154)-(  158).  It  i'ollow.s  from  the  discussion  in  Appendix  G 
that  eqs,  ( 154 ) -( l.'iS )  characterize  an  object  matched  in 
angular  extent  to  the  Instantaneous  f ield-of-vlew  of  the  scanning 
sensor,  i.e,,  a  "turget-like”  object.  As  illustrated  In  this 
figure,  the  point-source  approximation  nipg  provides  an  excellent 
approximation  to  for  target-like  objects. 
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NORMALIZED  THRESHOLD 


Figure  .19 

Similar  to  Fig.  .18.  except  that  the  field-ot’-view  is 
now  .scanned  across  a  "clutter- 1  ike”  object,  five  times 
larger  in  angular  e.xtent  than  the  "target-like"  object  of 
Fig.  18.  The  solid  curve,  dotted  curve,  and  dashed  curve, 
correspond  respectively  to  the  quantities  .m  ^  and  , 

As  illustrated  In  this  figure,  the  asymptotic  approximation 
!t\ .  provides  an  excellent  approximation  to  m,  for  clutter- 

t\  i) 

like  ob.jects. 
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Figure  20 

A  simple  adaptive  threshold  scheme  is  illustrated. 
The  transfer  function  exp(-27rf  T^)  introduces  a  delay 
of  seconds ,  ensuring  decorrelation  of  the  random 
processes  Y(t)  and  Y^(t)  (cf.  Fig.  21),  The  significance 
of  the  random  currents  is  seen  by  inspection  of  Fig.  1. 
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A  HiinpK'  iulapt;ivo  ihri'shold  .'U'ln'tiu*  1h  i  1  1  ust  rii  l  ml .  'I'lu' 
doLay  olonund.i.i  I'iioli  Lnt.rv?duo»‘  a  iU?lay  of  sih'oiuIh,  Tho 
f  loanoo  o  r  t,li('  random  v’urront.s  X(  t.  3 .  Y(.  t  )  .  and 
iH  soon  by  insjH'otlon  of  Flu'.  1.  Tho  blook  vlliiKrnm  shiuvn 
liori'  Is  aot.ually  Just  ono  part  of  tin'  tlirosliold  lU'oross  1  np' 
roiM'lvor  slunvn  tn  IMk-  I- 
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X(t) 


Finuro  23 

A  candidatu  adaptive  thr«^shold  scheme  is  illustrated  that 
tjeneralizes  the  structure  of  Fig.  22 .  The  threshold-establi.shiug 
approach  .shown  hero  is  realiaod  in  terms  of  a  tapped  delay  line 
with  2N  tup.s.  Once  again,  the  .significance  of  X(t),  Y(t),  and 
Y^(t),  follow.s  from  Fig.  1. 
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Figure  24 

The  tapped  delay  line  adaptive  threshold  scheme  shown 
here  generalizes  the  structure  of  Pig.  23.  The  transfer 
functions  H^^Cf),  n*l,2,...2.N  at  each  of  the  2N  taps  of  the 
tapped  delay  line  are  de-  i  variables,  chosen  to  mariimize 
the  receivers'  performance  against  a  particular  background 
scene,  or  set  of  background  scenes.  There  is  no  a  priori 
reason  why  the  various  delays  should  be  chosen  a.s  equal 
to  one  another,  other  than  for  fabricational  simplicity. 

More  generally,  additional  degrees  of  freedom  aro  incorporated 
by  allowing  these  delays  to  take  on  distinct  values. 
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APPENDIX  A  •“  Calculating  the  Average  Photocurront  from  BaokKround 
bat  a 


1 .  Objective 

The  objective  of  this  Appendix  is  to  show  how  the 
average  value 


m^it)  -  E{X(t)} 


of  the  random  current  X(t)  (cf.Fig.  1)  can  be  derived 
in  terms  of 

a)  high  spatial  resolution  radiometric  imagery  of 

the  infrared  scene,  in  the  optical  passband  of  interest 

b)  the  postulated  Modulation  Transfer  Function  (MTF) 
of  the  optical  train  of  the  model  system, 

c)  the  quantum  efficiency  and  physical  dimensions 
of  the  focal  plane  detectors,  and 

d)  the  focal  plane  scan  velocity. 

The  basic  result  to  be  derived  is  eq .  (A-14). 

2 .  The  Focal  Plane  Irradiance 

In  order  to  calculate  ra^(t),  it  is  first  necessary  to 
obtain  an  expression  for  the  average  irradiance  mg(x,y)  that 
would  be  measured  in  the  focal  plane  of  the  model  system 
in  terms  of  the  average  aperture  irradiance  mjjQ(.x,y).^ 


^The  aperture  irradiance  is  a  random  process  due  to  photon 
noise.  Thus,  following  the  usual  convention  of  random 
process  theory,  the  complete  ensemble  of  aperture  irradiance 
functions  is  denoted  by  a  capital  letter,  E(x,y).  A  single 
member  of  the  ensemble  is  denoted  by  a  lower  case  letter, 
o(.x,y).  The  capital  letter  E  is  used  to  denote  both 
irradiance  and  the  process  of  statistical  expectation;  however, 
the  context  of  usage  should  eliminate  any  possibility  of  con¬ 
fusing  which  of  these  two  meanings  Is  meant. 
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Tho  mean  aperture  irradiance  is  cluiraoteristic  oi  a 

particular  infrared  scene,  and  may  be '-■stlmatcd  liy  rnc.inH 
of  a  radiometric  Thermal  Iimn-fiUK  System^  (TiS)  of 

hij'her  spatial  resolution  than  tho  model  systtmi.  It  is  also 
highly  desirable  that  the  dwell  time  of  the  TIS  be  much 
larger  than  that  of  the  model  system,  since  the  analysis 
requires  knowledge  of  the  mean  irradiance  established 
by  averaging  over  tho  photon  fluctuation  statistics  of  thi' 
incident  light. ^ 

Tho  spectral  filter  chosen  for  use  with  tho  TIS  shi)uld 

match  tho  combined  TIS  optical  train/photodott'ctor  spectral 

response  to  that  of  tho  modol  .system.  Tills  i.s  nocessary 

bocuuse  there  is  no  way  to  re‘liabiy  calculate  tho  irradiance 

of  a  scene  measured  in  a  waveband  ' ^ j  ih  terms  of  tlu' 
irradiance  of  tlu'  sumo  soinio  mi'usured  in  a  dU'ft'rt'iu  spectral 


bund  A\., . 


Once  m,,^(x,y)  has  been  .specified,  the  "transformed  mean 


aperture  irradianee"  ^-‘o^  *  v ^  obtainiui  as 

(?  (  f  .  f  )  =  r  .y) ) 

X  y  EO 


I  A-  n 


■^A  TIS  with  a  dwell  time  mutehod  to  tho  dwell  time  of  tho 
model  sy.stem  would  actually  be  measuring  a  single  member 
e^(x,y)  of  the  ensemble  K  (x.y).  rather  than  the  expeeted 
vhlue  m„Q(x,y)  iif  the  apei^tiire  Irradianee,  Averaging  out 
the  photon  noise  by  either  increasing  the  TIS  dwi'll  time 
or  by  performing  friuac*  addition  decreases  tlu?  noisiness  of 
tho  TIS  image,  and  improves  the  goodness  of  the  TIS 
imagery  as  an  estimate  of  nipt  x.y). 


10') 


whero 


i  .  1  denotoii  the  Fourier  transform  of  the  braoki^ted 
quant  tty 

(f  ,f  )  are  spatial  frequencies  correspond  inn  to  the 
'  ^  space  coordinates  x  and  y,  respectively 

Employ ina  vector  notation  eq.  (A-1)  may  be  written  as 

The  transformed  focal  plane  Irradiance  (?( p  is  obtained 
by  mult  ip  I  y  inw  ^)by  the  Modulation  Transfer  Function  UlTF) 

that  characterizes  the  imane  blurring  offiH't  of  the  model 
systems'  optical  train: 

-  mtfq;)  (A-121 


whore  K  is  a  constant.  Assuming  the  convent  tonal 
o 

MTF  normal  I zat Ion : 

MTF  -  1 

conservation  of  energy  rc'quiros  that; 

Whore ; 

■  the  transmittance  of  tile  optics 
»  the  focal  lenqth  ratio  of  the  optics 


no 


Finally,  thi!?  mean  I’lioul  plane  Irrudiance  is  obtained  by 
Lakinj;  the  inverse  Fourier  transform  of^''(,^): 

iiijj  t  ^ ,|’ )  }  ( A-:i ) 

h  .  The  Timo-Vary  iutt  Average  Photi^ourrent 

The  response  itf  a  photovoltaic  detector  centered  at 
position  in  tlu'  ft)oal  plane  may  be  writtcMi  as; 


'^det 


(A-4) 


wliere  the  .surface  Integral  extends  ovt'r  iht'  ari‘a 

I 

of  the  detector,  and  is  an  clement  of  surface 
art'a  (cf.  Fl^'  A.l). 

The  focal  plane  Irr  ad  lance  and  the  moan  current 
in  cq,  (A-4)  are  onst'mble  averaf^c  values  over  the  photon 
f  1  uc  Lua  t  ioa  statistics  of  the  lip;ht.  For  the  purpose  of 
analysis,  the  electrical  freuuency  character ist Ic  of  the 
detei'tor  has  been  lumped  toKctnor  with  the  transfer  function 
of  the  pi\st-detector  filter.  Tlius  .  by  def  in  i  t  ion  .the 
current  respi>ns  1  vi  ty  K  ^  in  eq  .  (A-4)  is  Independt’n  t  of 
the  frequency  of  any  temporal  modulation  that  may  bo  ImposiHl 
on  the  trrudtancc. 


Ill 


A  "pupil  fuaotlon"  ) 

la  now  defined 

for  t  lu)  deto>v:tor 

such  that 

1 

>  V. 

do  t 

( A-fj  1 

1 

«  r' 

••1  "^dt.'t 

Aiisumii\v?  that  tUe  has  a  of  syuwutry,  so 

chat 

P(\*  “'I  -  Pf-V. 

and  that  tho  i’t>sponslvl  ly  Rj  Is  uni  form  on  tin'  ilotootor 
surface,  on,  tA— n  tioeoinos 

;U  ’"K^Von*  ' 

"'n^^o^  "  ;v.. 

where  A  is  the  Inflttlt  o  t\H'ul  plane.  Wlti»  the  following 
ohanne  of  variable 

-  V,,  V.' 

on,  (A-il^  beoonu's 


nt,.  1  r  ) 
\  \.o 
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(A-7) 


DefUUuK  th,»  t’ol  lowing  Lwo-di,nonHlonal  Fouriyr 


trun  t’oniiH  ; 


6'(4;)  -  K|in^(^)l 


C  A-H  ) 


I  -  F(P(^')} 

Eq  (A-7}  may  be  fcoiisf  with  the  aid  of  the  Fourier  e 
t  Ion  theort-in  ii8 


(  A-P ) 


onvol u- 


-  Rj  eil)} 

Mor«\)ver.  it  t'oUowH  from  oq.  iA-2)  that 


lA-lO) 


&'({)  -  MTFC,()  .  K 


(A-U  ) 


whia'o 


sPiiMuI  Fourl.'i'  trims  for 


orm  of  the 


nporniro  irnidianco.  and 


MTF(|^)  i.s  the  modulation  transfer  fun.' 


or  function 


of  the  optics. 


From  (A-IO^  and  (A- 1.1) 


CA-12  ) 


whore,  by  delMnition. 


MTF  ip  ;;  pp-)  MTFtp  •  K 


(A-i:n 
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For  tnttnitoaml  "point  detectors"  is  effectively  unity, 

so  that 

MTF '  -  MTF  . 

for  this  special  case.  More  generally,  the  effect  of 
having  a  finite  detector  size  can  be  accounted  for  by 
multiplying  the  MTF  of  the  optical  system  by  the  transformed 
pupil  function  of  the  detector,  as  in  Eq.  (A-13). 

Finally,  scanning  effects  are  included  by  allowing  the 
vector  in  Fig.  A.  1  to  trace  the  appropriate  trajectory. 

For  a  uniform  scanninL:  rate 

Ko  " 

where  t  is  the  time,  and  ^  is  the  focal  plane  scan  velocity. 

Eq,  (A- 12)  then  becomes 

m^(t)  -  Rj  /  MTF’C;^)  oxp  (j2TT^.Yt)d|  (A-IA) 

The  current  responsivlty  may  be  written  as: 

Rj  =  ne/hv  (A-15) 

wh  or  e 


h  “  Planck's  constant ( joule-soc . /photon ) 
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I 

I 

\ 


[Jtni  1  o-Mor  ■ _ j 

plui l  on  J 

'  h-j 


U  «  IMaiu 
V  ="  >'!'t  U'a  I  fr'.Hiut'nry 


n  =•  luuumiiti  o  t' t' lo  1  mu 


r  pliotoo  U'o  r  rooH  1 

'V  liu*  1  ilont  p’luuotiH  / 


®  t' I I'r  t. non  l  o  I'lini’Ko 


[noil  1  ombsl 
(» Uu‘  I.  noiij 


Hu'  liu-al  piano  i-sinin  V(»loolty  in  tni .  (,A->L4)  Is  ^'Ivon  by; 

I A-lUi 


\v 


VI  - 


WhiM'i' 


w  =«  t  lu'  do  t  oo  tail'  width  o 


“  syst  om  dwoll  f.  iitu'  soo 


It  may  ho  shown  t’luim  t'q .  (A- 10)  that: 

IV,I  -  2.  f-D  /T,.. 


t A~17^ 


whol’o  : 


^  t  lu’  I'lvul  lonuth  ratio  of  t.lU'  optics 


ditnons  1 1.1111  oss 


'\ip  “  ' ^iIi‘mt‘tor  of  tlio  optloal 
Tj,,  =*  syst  i'm  I'ramo  tlmo 


apor furo 


oni 
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. 

_ _ J 

The  focal  plane  irradiance  stationary  in  the 

coordinate  system  with  origin  0.  Vector  locates  the 
center  of  a  detector  of  area  For  scanning  sensors 

is  a  function  of  time. 
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APPENDIX  B.  NOISE  CURRENT  CORRELATION  FUNCTIONS 


The  object  of  this  Appendix  is  to  verify  eqs .  (73)-(76) 
which  are  needed  to  evaluate  eq.  (70)  for  the  threshold  crossing 
rate  ihj(  t ) . 

The  starting  point  for  this  discussion  is  the  linear  system 
input/output  relation  between  the  random  processes  X(t)  and 
Y(t)  (cf.  Figure  1): 

00 

Y(t)  -  /  du  X(li)h(t-u)  '  (B-l) 

—  TO 

It  follows  directly  from  (B-1)  that; 

00 

«  //  rtudX  h(t^-X)  hCtg-U)  (B-2) 

**  <*oo 

where  the  covariances  Cyy  s^hd  are  defined  by 

Cyy^^1'^2^  5  E{CY(tj^)-ray(t^)J  i:Y(t2)-my(t2)3  }  (B-3) 

CjjjjCX,).)  3  F,{CX(X)-m^U)J  CX(u)-mj^(y)J}  (B-4) 

and  where 

my(t)  =  E{Y(t)}  (B-5) 

m^(t)  s  E{X(t)}  (B-6) 

7i 
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With  the  definition 


Z(t)  =  Y(t)  (B-7) 

it  follows  from  eq.  (B-3)  that 

=YZ<*1’‘2>  ■ 

and 

'^Z2<'r'2)  ■ 

From  eqs.  (B-2),  (B-8)  and  (B-9) 

go 

,t^)  -  //  dvidA  h(t^-A)  3^  hCtg-U)  (B-10) 

•mOO  2 

00 

C22(ti,t2)  -  //  dudX  C^^CX.P)  3t  3^  hCtg-y)  (B-ll) 

-oo  1  2 

Setting  tj^=t2“t  in  eqs.  (B-10)  and  (B-ll),  and  noting  that 

C22(t,t)  =  o^ct)  «  EfZ^Ct)}  -  m^(t)  (B-12) 

it  follows  that 

00 

CygCt.t)  =  //  dydX  3^h(t-A)  (B-13) 

a^Ct)  =  //  dudA  a^h(t-X)  D^h(t-u)  (B-14) 
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An  expression  for  C  (A,u)  is  now  required  before  the  analysis 

XX 

can  be  carried  any  further. 

Setting  in  aq.  (B-2),  and  noting  that 

CyyCt.t)  -  a^(t)  -  E{Y^(t)}  -  m2(t)  (B-15) 

it  follows  that 

00 

al  (t)  -  //  dydA  a^(A.ia)h(t-A)h(t-u)  (B-16) 

y  _oa 

However,  an  adaptation  of  eq.  (4.3.13)  on  p.  115  of  ref.  (12) 
leads  to : 

00 

a^(t)  -  e  /  dy  h^(t--]X)  m(u)  (B-17)  • 

where  e  is  the  electronic  charge.  Consistency  between  eqs . 

(B-16)  and  (B-17;  requires  that 

Cj^j^(A,y)  -  e  nij^CM)  6(A-u)  (B-18) 

where  6(.)  is  the  Dirac  delta  function.  Covariance  functions 
like  eq.  (B-18)  are  characteristic  of  non-stationary  white 
noise. From  eqs.  (B-13),  (B-14),  and  (D-18), 

00 

CY2(t,t)  =  e  /  dy  h(t-y)  9^h(t-y)  tnj^(y)  (B-19) 
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2 


(B-20) 


a2(t)  -  e  /  du 

— oo  * 


m,(u) 


Noting  that 

h(t-ii)  3^h(t-u)  ■  i  3^  h^(t-u) 

it  follows  from  eqs.  (B-19)  and  (B-21)  that 

00 

CyzCt.t)  -  4  e  3^{  /  du  h^Ct-y)  m^(y)} 

mQO 

From  (B-17)  and  (B-22) 

CYz(t,t)  -  4  3^{a^(t)}  -  OyCt)  efy(t) 


Defining  r(t)  as 

r(t)  -  CygCt.t)  r'^y(t)o^(t)3 

it  follows  from  (B-23)  and  (B-2'4)  that 
r(t)  -  {*y(t)/a^(t)} 

Taking  the  expected  value  of  both  sides  of  eq.  (B- 
to  the  result 

CO 

ra  (t)  -  /  du  m^(y)  h(t-u) 


(B-21) 

(B-22) 

(B-23) 

(B-24) 

(B-25) 

leads 

(B-26) 
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Taking  the  time  derivative  of  (B-1): 

00 

y(t)  “  Z(t)  -  /  dy  x(vi)  3  h(t-vi)  (B-27) 

.00 

Taking  the  expected  value  of  both  sides  of  (B-27); 

00 

na^jCt)  -  f  du.tn^(u)  3^h(t-y) 

vOO 

1  >6.  I 

00 

m  (t)  -  3.{  I  du  rn  (y)  h(t-y)}  (B-28) 

^  ^  -w  ^ 

From  eqs.  (B-28)  and  (B-28) 

m^Ct)  -  rliy(t)  (B-29) 

Defining  the  convolution  operator  as  in  eq.  (77),  eqs.  (B-26), 
(B-17),  (B-20),  (B-29),  and  (B-25)  may  be  written  as: 


my(t)  - 

h(t)  •  m^^Ct) 

(B-30) 

a2(t)  - 

e  h^(t)  • 

(B-31) 

olit)  - 

ei;ri<t)}^  •  ra^(t) 

(B-32) 

in^(t)  - 

(B-33) 

r(t)  - 

i'cfy(t)/a^(t)} 

(B-34) 

Equations  (B-30)-(B-34)  are  .Identical  to  eqs.  (73)-(76). 
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Appendix  C.  Relationships  Between  FAR,  P^,  and  rHj 

1.  The  Relationship  Between  FAR  and  rUj 

The  coraplete  description  of  an  IRST  sensors’  perfonnance 
under  a  given  set  of  operational  conditions  requires  the 
simultaneous  specification  of  both  the  False  Alarm  Rate  (FAR)  and 
the  Probability  of  Detection  (Pjj)  for  a  "target"  within  the 
sensors'  field  of  view.  However,  both  "false  alarms"  and  target 
detections  are  manifested  as  threshold  crossings  by  the  signal 
processor.  Thus,  the  object  of  this  Appendix  is  to  relate  the 
traditional  IRST  performance  measures,  and  FAR,  to  the 
expected  number  of  threshold  crossings  mj  over  prescribed 
intervals  of  time. 

It  is  assumed  that  the  average  current  m^^Ct)  is  known  on 
an  interval  of  time  |t|  <  T/2 

The  expected  number  of  threshold  crossings  on  the  interval 
t|  <  T/2  is  defined  as  raj(0,T).  Defining  the  false  alarm  rate  as 
the  expected  number  of  threshold  crossings  per  "reference 
interval"  the  following  relationship  obtains  between 

FAR  and  mj  : 

far  -  /T)  mj(0,T)  (C-1) 


The  quantity  ra  (t)  may  be  specified  a  priori,  or  it  may  be 

A 

calculated  in  terms  of  the  radiance  of  a  particular  background 
scene  (as  discussed  in  Appendix  A). 
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For  example,  if  FAR  is  defined  as  the  average  number  of 
false  alarms  per  week,  is  set  equal  to  the  number  of 

seconds  in  one  week;  if  FAR  is  defined  as  the  average  number 
of  false  alarms  per  system  dwell  time,  then  is  set  equal 

to  the  dwell  time  (again  expressed  in  units  of  seconds). 

It  is  implicitly  assumed  in  eq.  (C-l).that  the  scene  under 
observation  does  not  Include  a  target,  so  that  each  threshold 
crossing  that  occurs  gives  rise  to  a  "false  alarm". 


2.  The  Relationship  Between  and  m,;  First-Order  Approxima¬ 
tion  ° 

Although  not  as  straight-forward  as  eq.  (C-l).a  relation¬ 
ship  between  and  can  also  be  established. 

AS  prelude  to  the  detleltlon  ol  P,,,  a  "deolslon  Interval" 

Tjj  is  first  defined.  The  interval  T^  is  presumed  to  bracket 
the  entire  period  of  time  during  which  the  current  Y(t)  manifests 
target-induced  fluctuations. 

Assuming  that  a  target  is  present  in  the  scene,  the  number 
of  threshold  crossings  that  occur  during  the  interval  is 
defined  as  the  integer  random  variable  J.  The  discrete 
probability  density  function  of  J  is  denoted  as  fj(J). 

The  probability  of  detection  P^^  is  now  defined  as  the 
probability  that  one  or  more  threshold  crossings  occur  during 
the  decision  interval: 

oo 

Pr,  -  ^  f/J)  (C-2) 

^  J-1 
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(The  likelihood  of  a  background-induced  crossing 
during  Tj^  has  been  neglected).  Unfortunately,  the  problem 
of  obtaining  a  formulation  for  fj  appears  to  be  quite 
difficult.  '  The  focus  of  this  paper  has  been  on  the 
development  of  formulations  for  the  expected  number  of 
threshold  crossings 


m-  =  E{J}  -  E  Jf,(J)  (C-3) 

j-1 


In  order  to  establish  a  relationship  between  and  mj, 
eqs .  (C-2)  and  (C-3)  are  written  as: 


00 

m,  -  fjd)  +  E  Jfj(J) 

Assuming  that  the  probability  of  two  or  more  threshold 
crossings  is  negligible  during  the  decision  interval  T^, 
eqs.  (C-4)  and  (C-5)  can  be  estimated  as: 

Pd  “ 


(C-6) 

(C-7) 
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It  follows  from  (C-6)  and  (C-7)  that 


Pjj  «  mj(0,T^)  (C-8) 

According  to  eq.  (C-8),  the  expected  number  of  threshold 
crossings  during  the  decision  interval  provides  a  good 
estimate  of  the  detection  probability  Pjj,  so  long  as  the 
probability  of  two  or  more  crossings  during  is  negligible. 

3.  The  Relationship  Between  and  m.:  Infinite  Contrast 
Point  Targets  ^  '' 

An  alternate  approach  is  now  taken  to  the  problem  of 
establishing  the  relationship  between  and  Pj^.  The 
discussion  that  follows  is  based  on  the  analysis  of  Section  IV, 
and  will  lead  to  the  same  result  as  before  -  namely,  that  the 
expected  number  of  threshold  crossings  during  the  decision 
interval  provides  a  reasonable  estimate  for  P^. 

It  is  shown  in  Section  IV  that  for  BLIP  sensors  and 
unresolved  targets: 

"■j  ■  ^ 

where  4>(.)  is  the  Gaussian  distribution  function 

♦  (u)  •“  (2tt)  *  /  exp  (-v^/2)dv  (C-10) 

«.oo 

and 

“min  -  “y<'mux> 
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where  t  is  the  time  at  which  the  mean  current  ™y(t) 
assumes  its  peak  value. 

Eq.  (C-9)  is  equal  to  the  exceedance  probability  of 
a  Gaussian  random  variable  having  a  mean  and  variance 
of  and  respectively.  As  such, 

eq.  (C-9)  is  a  relatively  straightforward  generalization 
of  the  procedure  conventionally  used  to  define  Pjj-  ^ 

Once  again,  one  is  encouraged  to  interpret  as  a 
reasonable  approximation  to 

4.  The  Relationship  Between  P^.  andm-:  Proposal  for  a 
Second-Order  Approximation^ 

The  first-order  approximation  to  discussed  in  Section 
0.2  Is  hs^sed  on  the  following  logic: 

a)  The  methods  discussed  in  Section  IV  of  this  paper  are 
used  to  obtain  a  value  for  mj,  the  first  moment  (i.e., 
the  average  value)  of  the  random  variable  J. 

b)  Following  eqa.  (C-6)  and  (C-7),  mj  is  used  as  an 
estimate  for  fj(l). 

c)  The  estimate  for  fj(l)  established  in  b)  is  substituted 

into  the  truncated  summation  for  Pj^,  eq.  (C-6). 

It  appears  only  reasonable  that  an  improved  estimate  could 

2 

be  obtained  for  P^  if  the  variance  Oj  of  J  were  known  in 

addition  to  the  mean  mj  of  J.  It  will  now  be  shown  how  knowledge 
2 

of  Oj  can  be  used  to  calculate  a  second-order  approximation  for 
Pjj  (compare  with  (C-6)); 
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Pd  ”  ^ 


(C-12) 


Unfortunately,  the  technique  described  in  this  section  for 

calculating  the  second-order  approximation  to  Pq  cannot  be  im- 

2 

plemented  until  a  formulation  for  a,  is  developed  analagous  to 

V  ■ 

the  formulation  for  mj  in  Section  III.  In  this  connection,  it 

is  noted  that  Bendat  has  derived  an  equation  for  the  crossing 

count  variance  of  stationary  processes.^  '  His  result 

(cf.  also  ref.  (23))  is  far  more  complicated  than  the  analagous 

eq.  (E-S)  for  mj.  Thus  a  generalization  of  Bendat 's  result  for 
2 

Cj  to  the  case  of  nonstationary  processes  and  stochastic  threshold 
functions  may  prove  to  b  )  a  difficult  problem. 

It  is  now  assunii'd  that  a  formulation  for  Cj  can  be  obtained, 
analagous  to  the  development  for  in  Section  III.  Analagous  to 
eq.  (C-3): 

00 

Oj  fj(j)|  -mJ  CC-13) 


Substituting  eq.  (C-3)  into  (C-13)  leads  to  the  following 


expression  for  Oj. 


-  fj(l) 


(C-14) 


+  4fj(2) 


l-fj(l)-fj(2)  +  Ej 

•  « 


where 


E,  =  E  rf 

j-3 


-2  [ 


f(l)  +  2f(2)]  E  Jf  -I  E  Jf| 
J  j-3  Lj-3  J 


(C-15) 
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Assuming  that  the  probability  of  three  or  more  threshold 
crossings  is  negligible  during  the  decision  interval, 
eqs.  (C-3)  and  (C-14)  are  approximated  as: 

(C-16) 

<C~17) 

Calculation  of  mj  according  to  the  method  of  section  III,  and 

2 

an  analagous  calculation  for  Oj,  enables  eqs.  (C-*ie)  and 
(C-17)  to  be  solved  for  approximations  to  fj(l)  and 
The  second-order  approximation  to  is  then  obtained  by 
means  of  eq.  (C-12). 

If  eq.  (C-12)  is  found  to  yield  an  appreciably  different 
result  than  eq.  (C-8),  third-order  or  even  higher  order 
approximations  to  may  be  required;  otherwise,  the  first- 
order  eq.  (C-8)  is  then  verified  as  a  good  approximation  for 


mj  «  fj(l)  +  2fj(2) 
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APPENDIX  D.  CRAMER  AND  LEADBETTER'S  EQUATION 
'  ^  14 1 

Cramer  and  Leadbetter'  •^have  derived  the  following 
equation  for  the  rate  of  positive  slope  zero  crossings  of 
a  non-stat ionary  Gaussian  random  process  Y(t)  (cf.  ref >(14), 
eq.  (13.3.2)): 

ih^^(t;0')  (l-r  )*  <t>\<^y/  {<J>(n)  n'l'(n)}  (D-1) 


where 


(D~2) 


and  the  quantities  m  ,  a  ,  a  ,  and  r  are  defined  as  in  eq.  (70). 

j  y  ^ 

On  p.  288  of  ref.  (14)  the  statement  is  made  that  the 
mean  rate  at  which  Y(t)  crosses  a  curve  yQ(t)  can  be  obtained 
from  eq.  (D-1)  simply  by  making  the  substitution 

"y  ' 

in  e<i,.'  Cd-'I). 

From  eqs.  (D-1)  -  (D-3) 

"(^)  4>(u){<t>(Pj^)  +  Pj^IKPj)}  (D-4) 


with  defined  by  eq.  (67). 

Eq.  (D-4)  is  similar,  but  not  identical,  to  eq.  (70).  The 
relationship  between  eq,  (D-4)  and  eq.  (70)  will  now  be  explained, 
and  it  will  be  shown  that  there  is  an  assumption  implicit  in 
(D-4)  that  limits  the  range  of  its  validity. 
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Two  alternative  formulations  for  the  curve  crossing  rate 
are  derived  in  Section  III.B. 

a)  Assuming  that  both  conditions  (49)  and  (50)  must 

be  satisfied  before  a  threshold  crossing  can  be  "counted' 
by  the  processor,  the  resulting  expression  for  the 
crossing  rate  is  given  by  eq.  (51). 

b)  Assuming  that  only  condition  (50)  need  be  satisfied 
before  a  threshold  crossing  can  be  counted  by  the  processor, 
the  expression  for  rhj  is  given  by  eq .  (53). 

It  is  now  proposed  that  a  third  type  of  processor  might 
be  constructed,  sensitive  only  to  those  threshold  crossings 
for  which 


y(t)  >  ^o(t)  (.D-6) 

I.e.,  it  is  now  assumed  that  the  processor  requires  the 
satisfaction  of  inequality  (49),  but  not  (50).  It  then 
follows  from  eqs.  (48),  (D-5),  and  (52),  that ; 


ihj(t)  -  !(«■)  -  KS^^) 


From  eqs.  (D-6),  (63),  and  (67)-(69): 


ihj(t 


)  -(-y)  (I- 


r^)^  <f)(u)  {4>(Pj)  +  Pj<KPj^)} 


(D~6) 


(D-7) 
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Noting  that  eqs.  (D-4)  and  (D-7)  are  Identical,  it  follows 

t 

that  assumption  (D-5)  must  be  implicit  in  Cramer  and 
Leadbetter's  prescription  for  calculating  the  curve  crossing 
rate . 

An  asymptotic  analysis  of  eq .  (D-4)  similar  to  that  of 
Section  III.D.3  leads  to: 

analagous  to  eq.  (115).  It  follows  from  (D-8)  that  the 

mean  crossing  count  m^^j^  is  sensitive  to  all  threshold 

crossings  for  which  rfi  >  .  For  example,  is  incremented 

y  o  oij 

by  unity  both  at  time  and  time  tg  in  Fig.  7.  Thus,  the 
asymptotic  character  of  reflects  the  requirement  (D-5)  , 
as  expected . 
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Appendix  E  -  Introduction  to  Scanning  Sensors  and 
Adaptive  Threshold  Detectors 

1.  Scanning  Sensors  Operating  Against  Uniform  Scenes 

a.  Current  Statistics 

When  the  sensor  of  Fig.  1  scans  across  a  uniform 
background  scene,  the  output  current  Y(t)  is  "statistically 
stationary".  The  meaning  of  statistical  stationarity  will 
now  be  discussed,  as  background  to  the  discussion  of  non- 
uniform  scenes  and  non-stationary  processes  that  follows 
in  the  next  section. 

It  is  assumed  that  the  sensor  is  scanned  and  re-scanned 
over  the  same  scene,  and  that  there  are  no  changes  in 
either  the  scene  or  the  sensor  from  one  scan  to  the  next. 

The  current  y<t)  during  the  course  of  any  one  particular 
scan  is  called  a  "stunple  function"  of  the  random  process  Y(t), 
(The  process  Y(t),  in  turn,  may  bo  thought  of  as  the  infinite 
ensemble  of  possible  sample  functions.)  The  current  sample 
function  obtained  on  the  n^'*’  scan  is  designated  y(t;n).  We 
now  .;onsider  a  particular  one  of  these  stunple  funct.ions,  y(t;l), 
depicted  in  Fig.  E.la. 


132 


Tho  tiiiu)  varlutionH  In  y(t:l)  hav«  thoir  origin  in  the  time- 
ot’-iirrivul  f  luo  taut  ions  of  tho  indlvtduai  photons  inoidont  on 
tho  dotoo  tor ,  *  Thus,  tho  f  luctiiatlons  in  yCt:i)  aro  indopwndont 
of  tho  scan  volooity  and  are  present  rogardloss  of  whether  tlio 
sensor  is  scantling  or  motionless. 


♦  If  Muoocsaivo  photons  arrived  at  uniform,  or  othi'rwiso 
predictable,  time  tnomni'nts  the  arrival  time  of  tho  "next" 
photon  would  ho  prodictablo  In  advance  -  in  violation  of 
tho  Hoisonborg  uncorta.inty  principle.  Thus,  even  nominally 
constant-intensity  sources  give  rise  to  i'and«.)m  photocurrentvS. 
This  kind  of  noise  is  often  called  "photon  fluctuation 
noise"  or  "quantum  noise".  When  the  quantum  noise  associated 
with  tho  background  light  is  the  dominant  noise  typo  in  the 
sensin' ,  the  seinsor  i.s  said  to  be  opi'rat  ing  in  the  "Ihickground 
Limited  Performance"  (BLIP)  reginu*. 


The  average  current  at  a  particular  instant  of  time  t ^  raay 
be  defined  as  the  "ensemble  average": 

1  N 

m  (t  )  -  lim  {IT  E  y(t  :n)l  (E-1) 

y  °  N+oa  n-1  ° 

In  order  for  Y(t)  to  be  a  stationary  process  it  is  necessary  that 
m  (t),  as  defined  in  eq.  (E-1),  be  independent  of  time.  Thus, 

if 

"■y(to)  ■ 

where  the  times  t^  and  tj^  are  totally  arbitrary  (c.f,  Fig,  E.l). 
Similarly,  the  mean- square  deviation  of  Y(t)  from  its  average 
value  (i.e.,  the  "variance"  of  Y)  may  be  defined  at  each  instant 
of  time  as 

o  1  N  r  T  „ 

c^(t)  =  lira  K  E  y<t;n)  -  m„(t)  (E-2) 

y  N-^oo  n-lL  ^  J 

2 

The  variance  a  ,  like  the  mean  m  ,  Is  independent  of  time  for 

y  y 

stationary  processes. 

Eqs.  (E-1)  and  (E-2)  are  satisfactory  for  illustrating  the 

concept  of  "ensemble  averaging";  however,  it  is  desireable  to 

have  a  different  means  for  actually  calculating  the  values  of 
2 

in  and  a  in  terms  of  standard  background  and  sensor  panuneters. 

y  if 


1.34 


It  La  first  neoesaary  to  define  the  averuno  value  of  the 
current  X(t)  (of.  IMg .  I): 


m 


X 


(E-3) 


where 


n  “  quantum  efficiency  of  the  detec 
e  ■  electronic  charge 


coulomb a 
electron 


elec 
ctor  plioT 


alectrona 

pTioTon 


Qq  average  background  photon  flux 


I 

X 


photona 


sec 


It  may  be  shown  that  tor  aonaors  operating  in  the  BLIP  regime 

the  mean  value  and  variance  of  Y(t),  originally  defined  hy 

eqa.  (E-1)  and  (E-2)  ,  can  be  calculated  in  terms  of  m^^  as  follows: 


-H(0)  (R-4) 

-  2emj^Af  (K-R) 

where  H(0)  is  the  zero-ordlnate  of  the  transfor  function  11(0 
(cf.  Fig.  1),  and  Af  is  t!u'  noise  haiulwldth  of  U(  f)  .  For  bandpass 
H(f),  H(0)"0.  It  follows  from  eq.  (D-4)  that  for  this  case. 

Since  the  scene  is  spatially  uniform,  the  average  photon 
flux  Qg  is  independent  of  time.  It  follows  from  eqs,  (li-h)  •• 

(E-6)  that  the  mean  m  and  variance  o  "  are  also  Independent 

w  y 

of  time,  Justifying  the  claim  of  stationarity  for  the  current 
Y(t). 
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Assuming ' that  the  transfer  function  H(f)  is  normalized 


as  follows: 


max  H( f)  -  1  (E-e) 

f 

the  noise  bandwidth  Af  In  eq.  (E~5)  may  be  calculated  from 
the  equation: 

Af  -  /  |H<f)|^  df  (E-7) 

o 

b.  Crossing  Rates  for  Constant  Threshold  Detection 
It  Is  now  assumed  that  the  signal  processor  of  Fig.  1  is 
implemented  such  that  the  threshold  y^  is  equal  to  a  constant, 

The  constant  threshold  y^^  is  depicted  on  the  sample  function 
plots  of  Fig.  E.l.  A  brief  outline  will  now  be  given  of  a 
method  for  calculating  the  average  number  of  times  that  the 
random  process  Y(t)  crosses  the  threshold  during  a  time  interval 
of  duration  T  seconds.  (The  relationship  between  the  mean  number 
of  crossings  mj  and  the  usual  search  set  performance  parameters 
and  FAR  is  discussed  in  Appendix  C.) 

The  e.xpected  number  of  crossings  mj(0,T)  may  be  written  in 
terms  of  the  crossing  rate  rti^  as: 

mj(0,T)  -  rfijT  (E-8) 

(  L3 ^ 

According  to  Rice.'  rfij  may  be  calculated  as 

rhj  =  exp  (-u^/2)  (E-9) 
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where 


with  Af  given  by  (E-7).  Also, 

u  =  (yo-%)/^  (E-11) 

with  lUy  and  given  by  eqs .  (E-4)  and  (E-5). 

The  quantity  u  defined  by  eq.  (E-11)  may  be  thought  of 
as  a  normalized  threshold  level. 

It  follows  from  eq.  (E-9)  that  the  expected  number  of 
threshold  crossings  drops  off  rapidly  as  the  threshold  level 
is  increased. 

2.  Scanning  Sensors  Operating  Against  Non-Uniform  Scenes 
a,  Current  Statistics 

It  is  now  assumed  that  the  sensor  of  Fig.  1  is  scanned  a 

number  of  times  over  the  same  non-uniform  scene,  and  that  there 

are  no  changes  in  either  the  scene  or  the  sensor  from  one  scan 

to  the  next.  A  number  of  sample  functions  of  the  resulting 

current  process  Y(t)  are  depicted  in  Fig.  15.2. 

Once  again,  the  ensemble  average  mean  and  variance  of 

Y(t)  are  defined  by  eqs.  (E-1)  and  (E-2).  However,  as  will 

now  be  discussed,  Y(t)  is  now  a  non-stat ionary  process,  i.e,, 

2 

m  and  a  are  functions  of  time, 
y  y 
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It  is  assumed  that  the  infrared  scene  encompasses  ref^lon 
of  blue  sky  and  clouds  where 

Qjjg  -  average  photon  flux  incident  on  the  detector  when 
blue  sky  is  being  observed. 

■  average  photon  flux  Incident  on  the  detector  when 
cloud  is  being  observed. 


The  photon  flux  Q  in  eq,  (E-3)  is  seen  to  be  a  function 


of  time:  Q  takes  on  the  value  when  a  clovid  is  in  the  field  of 
view,  and  a  different  value  Qgg  when  the  scanning  field  of  view 
includes  only  blue  sky.  Thus,  the  process  Y(t)  is  non-stationary 
when  the  scene  is  non-uniform,  since  the  mean  and  variance  of 
Y(t)  are  seen  from  eqs.  (E-3)  -  (E-8)  to  be  functions  of  time.* 
The  time-varying  mean  value  m^(t)  is  superposed  as  a 
dashed  curve  on  each  of  the  sample  functions  y(t)  depicted  In 
Fig.  E.2. 


♦For  the  present,  it  suffices  to  say  that  the  forms  of 

eqs.  (E-4)  and  (E-5)  indicate  that  a  time-varying  m^^  must 

2 

give  rise  to  time  varying  ra  and  o  .  However,  it  should 

V  v' 

be  noted  that  eqs.  (E-4)  and  (E-5)  are  only  strictly  valid 
for  stationary  processes,  i.e.,  for  time-invariant  m^^ .  Thus, 
for  a  given  time-varying  m  (t),  the  correct  values  of  m  and 

X  y 

2 

0  cannot  generally  be  calculated  from  (E-4)  and  (E-5), 

Jr 

Generalizations  of  (E-4)  and  (E-5)  strictly  valid  for  both 
stationary  and  non-stationary  processes  are  given  by  eqs.((i) 
and  (7). 
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b.  Crossing  Rates  Tor  Constant  Threshold  Detection 
The  performance  of  a  constant  threshold  signal  processor 


(  cf.  Fig.  1)  against  a  non-uniform  scene  can  be  characterized 
in  terras  of  the  quantity;* 


ffljCO.T)  -  /  rtij(t)dt  (E- 

where  T^  is  the  time  interval  jt]  <  T/2,  and 

where  mj(0,T)  is  the  number  of  times  that  the  current  y(t) 

can  be  expected  to  cross  the  threshold  level  during  the 


time  interval  T  . 

o 


Eq.  (E-12)  is  a  straightforward 
generalization  of  eq.  (E-8),  to  allow  for  the  possibility  of 
time-variable  threshold  crossing  rates  ihj, 

As  long  as  the  time  variation  of  m„(t)  is  slow  compared 

X 

to  the  time  variation  of  the  impulse  response  h(t)  of  the 
post-detector  filter  (of  .Fig.  1),  a  good  estimate  for  it>T(t) 


can  be  obtained  from  eq.  (E-9). 


*Cf .  Appendix  C  for  a  discussion  of  the  relationship  of 
raj(0,T)  to  the  usual  IRST  performance  parameters  and  FAR. 


12) 
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The  following  steps  are  then  followed  in  calculating  mj(0,T); 

(i)  The  tirae-varying  mean  current  is  derived  from  the 

time-varying  photon  irradiance  Q(t)  by  means  of  eq.  (E-3). 

(A  detailed  derivation  of  terras  of  the  background 

radiance  distribution  is  provided  in  Appendix  A. ) 

2 

(ii)  Estimates  of  ra  (t)  and  o  (t)  are  obtained  from 

X  y 

eqs.  (E-4)  and  (E-5).  (More  rigorously,  eqs,  (6)  and 

2 

(7)  may  be  used  to  obtain  ray(t)  and  (t).) 

(iii)  Eq  .  (E-ll)  is  evaluated  for  u(t). 

(iv)  Eq.  (E-9)  is  evaluated  for  rh,(t), 

(v)  Eq.  (E-12)  is  evaluated  for  mj(0,T). 

Numerical  examples  following  the  above  prescription 
typically  show  that  the  crossing  rate  function  il>j(t)  is  extremely 
sharply  peaked  <cf.  Fig.  E.3).  Consequently,  appreciable 
contributions  to  raj.(0,T)  only  accrue  in  the  near  neighborhood 
of  points  such  as  t„  in  Fig.  E.3.  It  is  shown  in  Section  III.D 

p  Cl 

that  eq.  (E-12)  may  be  approximated  as: 

mj(0,T)  =  Hp  (K-i:5) 

with  th.(t„)  obtained  from  eq.  (E-9).  The  quantity  6t„  is  the 

up  p 

effective  interval  of  time  during  which  m  (t)  remains  in  the  near 
neighborhood  of  its  peak  value,  from  the  standpoint  of  crossing 
rate  calculations.  An  expression  for  vSt  is  derived  in  Section 

J. 

III.D  (cf.  eq .  16  or  eq.  96). 
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The  Implications  of  eq.  (E-13)  for  system  performance 
are  illustrated  with  the  aid  of  Fig.  E.4. 

The  large,  relatively  slowly  varying  maximum  centered 
at  t  in  Fig.  E.4  is  presumed  to  have  its  origin  in  the 

H 

background  scene.  The  narrower,  lower  amplitude  spike 

centered  at  in  Fig.  E.4  is  presumed  to  be  due  to  a  "target". 

It  follows  from  eqs.  (E-13)  and  (E-9)  that  the  likelihood 

of  a  clutter-induced  threshold  crossing  grows  rapidly  os  tiie  threshold 

level  y^  in  Fig.  E.4  is  lowered.  It  follows  from  the  analysis 

in  Section  III, D  that  a  clutter-induced  threshold  crossing  (i.e., 

a  "false  alarm")  becomes  a  virtual  certainty  when  the  threshold 

level  actually  intercepts  the  mean  current  m  (t).  There  is 

¥ 

apparently  no  way  for  the  constant  threshold  processor  to 
detect  the  target  peak  at  t^  without  also  incurring  a  false 
alarm  arising  from  the  clutter  peak  centered  at  tp. 


C .  Crossing  Rates  for  Ideal  (CFAR)  Adaptive  Threshold  Detection 

The  performance  of  an  adaptive  threshold  processor  is 

illustrated  with  the  aid  of  Fig.  E.5. 

The  processor  is  presumed  to  have  some  means  for  deriving 

high-confidence  estimates  for  m  (t)  and  a  (t), defined  asm  sind  o  , 

y  y  y  y 

respectively.  When  WyCt)  is  "slowly-varying”  the  processor 
establishes  y„(t)  as: 


y^Ct)  -  fhy(t)  4.  a^ct) 


(E-14) 


The  threshold  is  interpolated  through  periods  of  "rapidly 
varying"  m^Ct)  by  means  of  a  smoothing  filter,  (of.  ,  for 
example,  the  neighborhood  of  t^  in  Fig.  E,5).  The  adaptive 
threshold  constant  K^,p  in  eq.  (E-14)  is  a  design  parameter. 

From  eqs.  (E-11)  and  <E-14) 

u(t)  -  K^,j,  j^Oy(t)/Oy(t)j  +  jftyCt)  -  my(t)j/ay(t)  (E-15) 


When  the  estimation  errors  are  sufficiently  small 


Oy(t)  “ 


(E-16) 


ftyCt)  «  '^y<^) 


it  follows  from  (E-15)  that 


u(t)  =  K. 


(E-17) 
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Thus,  u(t)  is  rendered  time-invariant  by  the  adaptive 
threshold  processor  when  there  is  no  rapidly-varying  target 
contribution  to  m^Ct).  When  there  are  no  targets  in  the 
scene,  thr  mean  current  >ny(i)  is  assumed  to  be  slowly- 
varying,  and  the  expected  number  of  crossings  during  the  time 
interval  may  be  calculated  from  eqs.  (E-9) ,  (E-12),  and 

(E-17),  as: 

2 

raj(O.T)  -  /  rtijQ  .-T 
o 

Since  both  iti^^  and  K^,p  are  time-  nvariant,  eq.  (E-18)  becomes 

2 

mj(0,T)  -  exp  (-K^,p/2)  T  (E-19) 

Eq.  (E-19)  has  the  same  form  as  the  crossing  rate  expression 
for  uniform  scenes,  eq.  (E-8). 

The  kind  of  processor  just  described  has  been  called  a 
Constant  False  Alarm  Rate  (CFAR)  processor,  since  the  threshold 
crossing  rate  is  now  independent  of  time,  i.e. ,  a  crossing  is  no 
more  likely  to  occur  when  scanning  a  region  of  non-uniform  back¬ 
ground  than  when  scanning  a  region  of  uniform  background.  E.g.  , 
with  reference  to  Fig.  E.5,  the  crossing  rate  is  now  no  greater 

at  t„  than  at  any  other  time. 

P 
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It  shoiiUi  bo  boriio  In  mind  that:  ttin  (’FAR  proconsDr  In  a 

non-i*oal  i/,ubl.o  idoul;  I  t  1ms  boi'n  ussnint^d  that;  t;ho  proooHHDr 

Is  abl<'  to  t'stlmato  tho  nuunt  I  t,i,os  iti^(t;)  luul  ^’^It)  to  us  hl>ch 

a,  proclslon  as  doslrod.  tlonorally,  orrors  la  tlu‘  tustlmiitod 

valuos  for  in  and  o  uro  unavoUiublo,  giving  rlso  t.i)  approalablo 
y  )r 

t tmo-doptMuloutu*  Iti  tjq.  (K-IR)  for  u(t,). 


KlKuro  K.l 

Tlwoo  siunplo  t'\iuotluns  vil'  tlu'  ruiuiom  ourrf'nt  Y(i)  uro 

shown  us  I'uiu't.lons  of  t.lmo.  Tho  o\»rront  Yft.)  Is  f  ukon  ut  tin' 

output,  of  t  l\o  post -dtitootor  fit  tor  us  tlvo  .‘U'i\sor  Is  sounnod 

ov».u*  II  uniform  scoui'  (of.  Ktj;.  1).  Tlu'  snmplo  funotlons  nro 

dt's  I  ^^nuf  t‘d  y(t;:n),  ti- 1  ,  U  .  ;i ,  .  .  .  Tho  sivmplt>  fuiiot  U'us  hlspluy 

random  f,  l.mo  variations  oausotl  by  f.tmo-of-arrl  vul  fUiotuatlons 

of  tho  Individual  photons  Inoldont  luv  tho  dotootor.  Slnoo 

tho  ourront  Ylt)  l.s  stationary,  tho  onsomhlo  avoran'o  nu'iin 

valuo  (vurlanoo)  dofltu'd  by  oq .  (K-M  (oq,  K-f!1  Is  t  l\o 

.samo  at  t  Imo  t  a.s  at  tlmt*  t,,  whon.'  tlnu's  t  ^  and  t 
o  1  o 

arbitrary . 
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Three  sample  functions  of  the  random  current  V(t)  are 
shown  as  functions  of  time.  Thi.s  figure  is  similar  to  Fig. 

E.l,  except  the  sensor  is  now  presumably  scanned  over  a  non- 
uniform  scene.  The  time-varying  ensemble  average  m,^(t)  of 
Y(t)  is  shown  as  a  dashed  curve  superposed  on  each  of  the  three 
depicted  sample  functions  (solid  curves).  The  ensemble 
average  is  still  defined  by  eq.  (E-1);  however,  the  fact 
that  Is  now  a  function  of  time  Implies  that  Y(t)  is  now 
a  nonstationary  random  process. 
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my(t) 

Vo 

a) 

\ 


mj(t) 


b) 

Figure  E.3 

Part  a)  is  an  illustrative  plot  of  m  (t)  Vs.  t,  where 

y 

niy  is  the  mean  value  of  the  filtered  current  Y(t)  (cf,  Fig, 

13).  Also  shown  is  a  constant  threshold  current  ,  lying 

above  the  peak  value  of  m  .  The  function  m  (t)  takes  on 

y  y 

its  peak  value  at  the  time  t  . 

P 

Part  b)  is  a  plot  of  the  threshold  crossing  rate 
corresponding  to  the  threshold  y^  and  mean  current  m^Ct) 
of  part  a).  The  entirf)  contribution  to  the  crossing  rate 
integral,  eq.  (E-12),  accrues  in  the  very  near  neighborhood 
of  tp. 
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to 

Figure  E.4 

This  figure  Illustrates  a  critical  shortcoming  of 
constant  threshold  processing.  The  slowly  varying  maxima 
centered  at  tp  presumably  has  its  origin  in  the  non- 
uniform  background  scene.  The  narrower,  lower  amplitude 
spike  centered  at  t^  is  due  to  a  target.  The  likelihood 
of  a  false  alarm  (i.e.,  a  clutter-induced  threshold  crossing) 
grows  rapidly  as  the  threshold  level  y^  is  reduced.  There 
is  no  way  for  the  constant  threshold  processor  to  detect  the 
target  peak  at  t^  without  also  incurring  a  false  alarm 
arising  from  the  Clutter  peak  centered  at  t  . 

A  plot  of  the  threshold  crossing  rate  ihj(t)  corresponding 
to  this  figure  would  show  that  the  probability  of  a  false 
alarm  is  far  greater  at  time  t  than  akt  any  other  time. 
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Figure  E.5 


This  figure  illustrates  an  important  potential  advantage 

that  accrues  to  the  use  of  adaptive  threshold  processing. 

The  mean  current  m  (t)  is  the  same  as  for  Fig.  E.4,  with  a 

clutter  peak  centered  at  time  tp ,  and  a  lower-amplitude 

target  spike  centered  at  time  t^^ .  The  adaptive  threshold 

y  (t)  is  presumably  able  to  accurately  track  the  slowly 
o 

varying  background  signal,  but  not  the  more  rapidly  varying 

target  signal.  Thus,  target  detection  is  assured,  while 

the  probability  of  a  false  alarm  is  kept  acceptably  small. 

As  contrasted  with  the  situation  of  Fig.  E.4,  the  probability 

of  a  false  alarm  is  now  no  greater  in  the  neighborhood  of 

time  t  than  at  any  other  time. 
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APPENDIX  F.  THE  ETX'.EWOK'ni  tT«RE(.TION  TO  rii. 

J 

The  objective  of  this  Appendix  is  to  estimate  the  relative 
error  in  ifij  due  to  the  assumption  of  a  strictly  Gaussian 
probability  density  function  for  the  current  Y(t)  (of.  Fig.  1); 

-I 

t/y)  »  )  (’'-1' 

where  iK  • )  is  the  Gaussian  density  function,  m  is  the  mean  value 

y 

2 

of  Y(t),  and  G  is  the  variance  of  Y(t),  defined  respectively 

by  eqs.  (F-3),  (56),  and  (58). 

It  will  be  shown  that  the  assumption  of  non-stationary 

Gaussian  statist ic.s  for  the  filtered  current  Y(t)  is  totally 

Justifiable,  for  typical  sensor  parameters  and  background  radiances. 

The  starting  point  for  this  analysis  is  a  truncated 

f  13  ^ 

Edgeworth, series  expiwsion  ^  “^of  f^Iy),  as  appropriate  for 
describing  Poisson  shot  noise  processes  Y(t):* 

fy(y)  “  '^(u)  -  (Ay/GOy'S  (F-2) 


41 

strictly  speaking,  Y(t)  is  a  conditional  Poisson  shot  noise, 
with  u  wideband  thermal  level  density  function.  However,  it  Is 
well  known  that  such  processes  may  be  modeled  as  constant  level 
Poisson  shot  processes  when  the  response  time  of  the  detector 
is  very  Icng  compared  with  the  coherence  time  of  the  incident 
light,  (of.  Uef.  (12)) 
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whore 


•Ku)  “  (2n)“^  e.\p  (-u“/2) 

(F-;3) 

-  d“  4>(vi)/du'^ 

(F-^) 

u  ■  (y-my)/ay 

(F-5) 

niy  -  h(t)  w 

(F-6) 

-  e  h^(t)  mj^(t) 

(F-7) 

-  e“  h^(t)  «*  m^(t) 
y 

(F-8) 

The  quantities  e,  h(t),  and  in  eqs.  (F-6)-(F-8)  an?  defined 

foiiowinK  oq .  (76).  Eqs.  (F-7)  and  (F-8),  the  variance  and  third 
semi-invariant  cf  the  nonstut ionary  prooe.ss  Y(t).  are  adapted 
from  Ref.  (12). 

From  eqs.  (F-2)-(F-4) : 

fy(y)  -  >^(u)  (1  +  Ky)  (F~9) 

where  the  error  term  o ^  is  Riven  by 

-  (A  /0a  {u(u--:n} 


( F-LO) 


According  to  eq .  (F-10),  the  relative  error  increases 

3 

as  u  ;  however,  it  follows  from  eqs .  (F-3)  and  (F-9)  that  the 

density  function  f  Is  of  negligible  amplitude  for  large  values 

¥ 

of  u.  It  will  now  be  assumed  that  values  of  u  greater  than 
eight  are  of  no  importance,  since 

tAy)  -  cr{exp<-32)>  -  efdO"^'^),  u-8  (F-11) 

¥ 

In  order  to  establish  a  reasonable  upper  bound  for  the 
value  u«8  will  be  used  In  eq.  (F-IO),* 

Cy  <  500  (Xy/6c;y2) 

Also, 

^^y  max'^'^y  rain^ 


(F-12) 


(F-13) 


where,  from  eqs.  (F-7)  and  (F-8): 


a^j  •em  .  /  h^(t)dt 

y  min  x  rain  „  v 

¥  _  AQ 


(F-14) 


X  ■  ra  /  h^(t)dt 

y  max  x  max  vvyv*., 

>  M..CD 


(F-15) 


From  (F-14)  and  (F-15); 


max  °y  min 


^'"x  raax/"'.x  min^^"'x  min 


/e)“*  X 


(P-16) 
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where 


*  3  1  r 

“  2 

-3/2 

x 

s  /  h‘‘(t)dt 

/  h  (t)dt 
.00 

(F-17) 

The  number  of 

photoeventa  per 

■  “x  min/® 

second  Is 

now  defined  as ; 

(F-18) 

Assxunlng  that  the  scene  brightness  doesn't  change  by  more  than 
a  factor  of  two  from  one  dwell  time  to  another; 

max/  "‘x  min>  2  (F-19) 

From  eqs.  (F-12),  (F-13),  (F-16),  (F-18),  and  (F-19): 

e  <  200  v“^x  (F-20) 

¥ 

It  is  easy  to  show  that  for  a  rectangular-shaped  h(t) 
of  duration  eq.  (F-17)  becomes 

X  -  Tp"*  (F-21) 

where  tp  is  the  system  dwell  time.  From  (F-20)  and  (F-21).' 

<  200(VTp)"^  (F-22) 
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For  a  current  of  100  namps  and  a  dwell  time  of  160  nsec. 

(of.  oqa.  (154)  and  (G-14))th0  number  of  photoelectrona  per 
dwell  time  is  given  by: 

-7  -1 

10  ooul.aec 

^  ^  ^  ■1,1,  .«■ 

VTjj  “  1.6x10  ooul.  electron  (1.6x10  sec) 

I  .e . , 

vTj3  -  10®  (F-23) 

From  eqs.  (F-22)  and  (F-23): 

Ey  <  0.02  (F-24) 

In  general,  eq.  (r-21)  greatly  over-estimates  the  value 
of  x>  aod  the  corresponding  upper  bound  for  e  ia  much  smaller 
than  indicated  by  eq.  (F-24).  This  la  because  h(t)  la  the 
impulse  response  of  a  bandpass  filter,  taking  on  negative  values 
as  well  as  positive  values.  Thus,  there  is  cancellation  in  the 
numerator  of  eq.  (F-17). 

It  will  now  bo  shown  that  the  relative  error  in  fy(y) 
can  be  used  as  an  estimate  for  the  relative  error  in  tfij.  This 
point  is  most  easily  illustrated  by  employing  eq.  (D-6)  for  rh^. 
However,  either  of  the  crossing  rate  expressions  (Rl)  or  (63) 
could  be  used  instead,  at  the  expense  of  complicating  some  of 
the  equations, 
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From  eqs.  (D-6)  aad  (52): 

00 

«hj(t)  -  ;  (yQ.z)  dz  CF-25) 

For  the  numerical  examples  that  have  been  worked  thus  far 

(such  as  those  of  Section  V),  It  has  been  found  that  the 

* 

correlation  coefficient  of  Y(t)  and  Y(t)  la  tnuoh  smaller 
than  unity; 

r(t)  <<  1  (P-26) 

where  r(t)  is  defined  by  eq.  (B-24).  Since  it  has  Just  been 
shown  that  Y(t)  has  very  nearly  a  Gausaiar  density  function,  and 
since  uncorrelated  Gaussian  variates  are  statistically  Independent, 
it  follows  that 

^y^(yo,zi)  “  lyCy^)  t^(z)  (f-27) 

From  eqs.  (F-25)  and  CF-27) ; 

OO 

thj(t)  -  fy(yo)  /  f^(z)cVz  (F-28) 

^o 

Substituting  eq.  (F-9)  Into  eq .  (F-28),  it  is  oeen  that  e ^ 
provides  an  estimate  for  the  relative  error  in  thj,  as  well 
as  the  relative  error  in 
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APPENDIX  G.  PARAMETER  CHOICE  RATIONALE  FOR 
EQUATIONS  ('l49y  AND~ { 150) 

The  object  of  thi.s  Appendix  is  to  describe  a  rationale 
for  choosing  the  background  and  filter  parameters  ,  Xj^,  t, 
f^,  and  t^  (cf.  eqs.  (149)  and  (150)). 

It  is  common  IRST  design  practice'  '  to  choose  the 
noise  bandwidth  of  the  post-detector  filter  as  follows; 

Af  -  (2t.j.)"^  (G-1) 

where 


Af  *  noise  bandwidth  of  the  post-detector  filter 
Tip  ■  the  time  required  for  the  scanning  sensor 
to  sweep  ucro.ss  a  target . 

It  is  not  difficult  to  show  that 

Tp  «  tXpTp/2Tr  (G-2) 

where 

Up  ■  the  angular  extent  of  a  typical  target,  at  the 
desired  detection  range  (radians) 

Tj.  ■  the  system  "frame  time,"  i.e.,  the  time  required 
for  the  .sensor  to  complete  a  full  360*  azimuthal 
sweep . 
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It  may  be  shown  that  for  the  transfer  function  H(f)  of 
eq.  (151),  the  noise  bandwidth  is  given  by: 

Af  -  1^1  +  exp  (-ZTTf^^t^^^ 

Assuming  that 

,  H(0)  -  2  exp  (-TTfo^t^^)  «  1 

It  follows  from  eq.  (G-3)  that 
Af  «  (t^/2)"^ 


Eq.  (G~6)  is  the  Justification  for  referring  to  as  a 
bandwidth  parameter"  of  the  post-detector  filter. 

From  eqs .  (G-1)  and  (G-D); 

"  ■'t  ^ 

Assuming,  for  example,  that 

a,j,  *  1  mrad. 

Tp  ■  1  sec. 


it  follows  from  eqs.  (G-2)  and  (G-6)  that 


and 


T,j,  ■  159  psec. 
t^  "  225.  psec. 


(G-3) 

(G-4) 

(G-5) 

"noise 

(G-f.) 

(G-7) 

(G-8) 

(G-9) 
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Assuming  (somewhat  arbitrarily)  that 


H(0)  -  iO"^  (G~10) 

it  follows  from  eq.  (G-4),  (G-9) ,  and  (G-10)  that 

-  6.91  kHz  (G-10) 

Thus,  a  rationale  has  been  provided  tor  choosing  the  filter 
parametors  and  t^  of  eq.  (149). 

To  complete  the  problem  specification,  it  remains  to  choose 
values  for  the  parameters  of  eq.  (150): 


■;  «  tho  time  roquii’ed  for  the  sensor  to  scan 
ucros.s  the  object 

"  background  brightness  parameter  (amperes) 
Xj  "  object  contrast  panuneter  (dimensionless) 


Analiigou.s  to  eq.  (0-2): 


T  ■  a^^Tp/27r 

whtn'o 


(G-1.1) 


«  object  angular  e.'ttent  (radians) 
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The  form  of  eq .  (150)  is  general  enough  to  describe  both 
"target-like"  and  "clutter-like"  objects.  Target-like 
objects  are  defined  by: 


(G-12) 


Clutter-like  objects  are  defined  by 


a  ^  u™ 
o  T 


(G-13) 


It  follows  from  eqs.  CG-2),  (G-8),  and  (G-11)  that  the 
parameter  value 


T  ■  159  usee. 


(G-14) 


corresponds  to  a  target-like  object 

The  value  for  in  eq.  (150)  is  obtained  as  follows: 


X  «  R  P 

1  opt 


(G-15) 


where 


»  uniform  background  current  (amperes) 

Rj  «  current  responsivlty  ( lunperos/watt ) 

^opt  ”  optionl  power  incident  on  a  single  detector 
( wattvS) 


The  optical  power  in  eq .  (G-15)  is  calculated  as: 


^'opt  “  Si  ^lp 


(G-liO 


160 


where 


-2  ~i 

Lq  ■  in-band  background  radiance  (watts  cm  .  sr  ) 

2 

»  area  of  the  telescope  aperture  (cm  ) 
ua  ■  instantaneous  field-of-view  (sr) 

The  current  responsivity  Rj  in  eq.  (G-15)  is  calculated  as: 

Rj  -  Sv  (G-17) 

where 

n  “  quantum  efficiency  (photoelectrons  • 

-1 

(incident  photon)  ) 

e  -  electronic  charge  (coulombs  (electron)”^) 
h  “  Planck's  constant  (joule  sec  (photon)”^) 

V  »  optical  frequency  (sec.~^) 

It  follows  from  eq .  (G-17)  and  the  niunerical  values  of  e  and  h 
that : 

Rj  -  O.anA  (G-18) 

where 

X  =  average  optical  wavelength  (pm) 
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It  follows  from  eq .  (G-15),  (G-16),  and  (G-18)  that 


*  134  nomps.  (G-19) 

for  the  following  assiuned  parameter  values; 

2 

■  500  cm 
ap 

CD  “  a  *10  sr. 


and 


n  ■  0.0 

\  ■  4.0  um 

-2 

Lq  *  0 . 14  mw  cm  sr 


(G-20) 


The  value  of  background  radiance  given  by  eq.  (G-20)  was  obtained 
by  integrating  Pig.  1  of  Ref.  (28)  over  the  wavelength  band 
3~5  urn. 

The  assumed  value  of  peak  object  contrast,  eq .  (155),  was 
chosen  arbitrarily.  It  is  not  difficult  to  calculate  the  value 
of  apparent  contrast  radiant  intensity  corre.sponding  to  given 
values  of  and  Xj^. 
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